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ABSTRACT 


The results of a semi-analytical and experiment a], 
study of adiabatic tube flow and an analytical and experi- 
mental investigation of the thermal entrance region for gas 
flowing through electrically heated circular tubes are 
presented. Emphasis is placed on the low Reynolds number 
turbulent flow regime--def ined as fully turbulent flow at 
bulk Reynolds numbers from 3,000 to about 15,000. 

Adiabatic air velocity and friction data and localized heat 
transfer measurements for air and helium, at low heating 
rates, are presented for this range. 

The adiabatic data were obtained in a 1.6l inch ID 
tube for flow at bulk Reynolds numbers from 3 ,000 to 
15,000. A continuous, Reynolds number-dependent, profile 
is developed from the data by using a modification of 
Reichardts ' wall and middle law eddy diffusivity expres- 
sions . The velocity profile satisfies continuity. It is 
valid for all Reynolds numbers in excess of 3,000 for which 
the flow is fully turbulent and the Blasius friction factor 
expression is valid. 

The thermal entrance problem for a fully developed 
velocity profile is solved analytically by the method of 
Sparrow, Hallman, and Siegel. The solution is based on the 
profile developed from the velocity study. Tabular values 

xi 
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of the eigenvalues and normalized Nusselt numbers for gases 
are presented for a range of Reynolds numbers from 3 >000 to 
50,000. The axial variation of Nusselt number is found to 
be correlated by 

N “ "p " 1 

- = 1 «• 0.8(1 + 70,000 Re ‘'Mjj) 

u oo 

to within per cent for x/D 2 * 2. The fully developed 
value agrees with the Dittus-Boel ter correlation, 

v r\ nm n 0.8 0.4 

Nu = 0.021 Re Pr 

oo 

2 

For the eigenvalues , ^ , and the associated constants, A • 

correlations of the form 


-b. -d, 

i » n + q p^. ^ 

1 ,n 


) ~ = A. Re + C, Re 

An 1 , n 


-b 


2 ,n 


-d. 


A = -A,, Re “’** + Re 

n 2 , n 2 , n 


in 


are obtained. The coefficients and powers are presented in 
tabular form. 

Heat transfer data are presented, primarily for 
helium, for the low Reynolds number turbulent range. A 
one-quarter inch, resistively heated, vertical, circular 
tube was used for the study. The data cover an axial range 
from 1.2 to 96 diameters; wall - co-bulk temperature ratios 
vary from 1 to 1 . 4 . In the low Reynolds number turbulent 
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regime, these data clearly support the present analytical 
solution rather than the prediction obtained by applying 
the eddy diffusivity distribution used by Sparrow, Mailman, 
and Siegel. 


CHAPTER J 


INTRODUCTION 

At present, the treatment of turbulent flow in 
circular tubes has been adequately developed for high 
Reynolds number flow only. When these treatments tare 
extended, the predicted velocity profiles (Jo not agree with 
data in the low Reynolds number range. Resulting friction 
factor predictions are compared to experimental correla- 
tions in Figure 1. Discrepancies exceed twenty per cent. 
Further, when heat transfer predictions for the downstream 
Nusselt number are based on these adiabatic velocity 
profiles, they overpredict the available data by as much 
as fifty per cent as the Reynolds number is reduced. 

In this work the term "low Reynolds number turbu- 
lent flow" will refer to fully turbulent flow for the tube 
Reynolds number range from 3000 to 30,000, i.e., the region 
of the discrepancies just mentioned. The criterion for 
’’fully turbulent flow" is an int ermit tency factor--the 
ratio of the time the velocity is fluctuating to the total 
time of measurement — equal to unity at the location of 
interest. An "asymptotic" or "universal" profile refers 
to a non-dimensional , Reynolds number-invariant profile 
which is valid for large Reynolds numbers. For the flow of 
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an incompressible fluid with constant viscosity, thermal 
conductivity and specific heat, it can be shown that 
ideally the velocity profile and normalized temperature 
profile approach invariant distributions as the flow 
progresses down a channel under constant wall heat flux. 
Fully' developed conditions are reached when the invariant 
axial profiles are approached. For turbulent flow, this 
condition usually may be considered to be attained within 
fifty diameters or less after the start of heating. There- 
after, changes in the friction and heat transfer coeffi- 
cients are indiscernible for practical purposes, and the 
fluid bulk temperature and the tube wall temperature then 
both increase linearly. The section of the channel wherein 
the temperature profile adjusts to the fully established 
value from the uniform value prior to heating is called the 
thermal entry region, or the thermal entrance. 

With a solution for a boundary condition of 
constant wall heat flux from the start of heating, problems 
involving any axial variation of wall heat flux can be 
handled by the method of superposition, provided that the 
governing energy equation remains linear. This equation 
will be linear in temperature in the limiting case of low 
heat transfer rates because fluid properties can then be 
considered constant. To date, no thermal entry solution is 
available for low Reynolds number turbulent flow. 
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Currently there is an increasing interest in the 
problem of rel aininariza tion of a gas which is being heated 
strongly. Such conditions may exist in solid core nuclear 
rockets and in the fuel ducts of engines designed for 
hypersonic flight. But no basis will exist for determina- 
tion of the difference between behavior with the strong 
fluid property variation, accompanying the high heating 
rates, and "normal" hehavior--un.l ess a thermal entry solu- 
tion is made available for the same Reynolds number range 
for the low heating rate conditions . 

To overcome the foregoing discrepancies, this study 
was undertaken to provide accurate predictions for (a) 
adiabatic tube flow emphasizing the low Reynolds number, 
turbulent flow regime, and (b) local heat transfer 
characteristics in the thermal entrance and downstream 
regions (based on the results of the flow study) . The flow- 
regime of interest is the bulk Reynolds number range from 
3000 to 30 , 000 , with the main emphasis placed on the 3000 
to 15,000 interval. The adiabatic velocity profile data 
presented cover the latter flow range and confirm earlier 
trends indicated by isolated data in this particular 
regime. Through inclusion of the Reynolds number as an 
additional parameter, an adiabatic velocity profile has 
been developed that adequately correlates the data in tnis 
regime, as well as for Reynolds numbers in excess of 
30 , 000 . 


; » ■ 
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Based on the improved description of the adiabatic 
velocity profile, the energy equation, a second order 
partial differential equation, was solved. r nder the 
idealizations of constant wall heat flux, fully developed 
velocity profile, and constant fluid properties, an 
eigenvalue thermal entrance solution was obtained. It 
yielded temperature profiles and local Nusseit numbers as 
functions of axial distance and bulk Reynolds number. 
Experimental data, for low and moderate heating rates for 
gas flow in smooth, circular, electrically heated tubes, are 
presented for comparison with the analysis. These data are 
extrapolated to zero heating rate results to correspond 
with the idealizations of the analysis. 

This study is expected to be directly applicable to 
the heat transfer problems of resistojets, nuclear powered 
space vehicles utilizing gaseous propellents, proposed 
central station and marine gas cooled nuclear reactors, and 
to heated duct flow, in general. It may find use in the 
treatment of gas cooled nuclear reactor start up, shut down 
and partial loss of flow accidents. The results can also 
be used in electrical heating applications and for circular 
tube heat exchanger analysis. 
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PREVIOUS API A HAT I C STUDIES 

An adequate representation of the velocity profile 
is a necessity in the treatment of internal heat transfer 
problems. As the first stop in the development of such a 
representation, the pertinent literature is reviewed. 
Emphasis is on the low Reynolds number turbulent range and 
deviations from asymptotic generalizations • 

For circular tubes, the shear stress can be 
expressed as 

T = +PuW T (2-1) 

Two essential ly r equivalent approaches have been used in 
trying to express the Reynolds stress, u’v’, in terms of 
mean fluid properties: the eddy r diffusivity approach, and 

the mixing length approach postulated by r Prandtl (2). The 
literature abounds with velocity profile formulations, 
based on innumerable variations of the two approaches, for 
the high Reynolds number, or asymptotic, regime. Semi- 
empirical. models have been employed for the eddy diffusivity. 
Integration of the momentum equation, containing the 
resulting diffusivity expressions with a constant or linear 
variation employed for the shear stress, yields a universal 
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velocity profile. The constants ol' the profile are deter- 
mined from experimental high Reynolds number u + versus y + 
velocity profile data. The unfamiliar reader is referred 
to the texts by Knudson and Katz (l), or Kays (2), or the 
critical literature survey on internal turbulent flow by 
Cess (3) for more complete details of existing treatments 
for high Reynolds number turbulent flow. 

For high Reynolds number flow, the experiment al 
data can be represented to within ten per cent by any of a 
number of u + versus y + velocity correlations, each inde- 
pendent of bulk Reynolds number. However , data by Sonecal 
(*t). Rothfus and Monrad (5)> end Page and his co-workers 

(6) show the normalized velocity profile is not invariant 
below Reynolds numbers of about 30 , 000 . 

To date, only r five studies have attempted to 
include the low Reynolds number dependency. Rothfus and 
Monrad ( 5 ) indicate that a unique correlation is obtained 
if u V^/u^ versus y u^/V^ is plotted rather than u versus 

y + . However, the lack of an adequate prediction for u,„ 

X 

reduces the usefulness of this correlation. Gill and Scher 

(7) i extending the reasoning of Van Driest (8), considered 


the mixing length as dependent upon the Reynolds number. 
Their eddy diffusivity expression was phrased to yield a 
velocity profile that varies smoothly' from a parabolic 
shape at Re = 1,800 to agreement with Nikuradse's data at 
Re = 100,000. Pai (10) approximated the velocity profile 



a 

with a second order polynomial, with coefficients also 
chosen to yield a step-wise continuous variation from the 
parabolic laminar profile to the asymptotic high Reynolds 
number universal profile. However, Rotta's (17) experimental 
observation of a continuous axial variation from laminar to 
fully developed turbulent flow, while the diameter Reynolds 
number remains constant during transition, invalidates any 
single curve extrapolation to the laminar profile as in the 
last two formulations. for an entrance region analysis 
Hanrntty and Johnk (ll) used Deissler's (12) high Reynolds 
number velocity expression for y* < 26 and a power law of 
the form Cy* in the turbulent core, with C and M ns 

functions of Reynolds number. The lowest Reynolds number 
employed in their invos tigation was 18 , 000 . 

McEligot, Orm*l;$ •* and Perkins (9) applied a two 
layer treatment with a variable laminar sublayer thickness 
as a basis for a downstream heat transfer analysis. Their 
continuous, Reynolds number dependent velocity profile was 
based on Senecal's data (*t) at Reynolds numbers of 3)000 to 
*1,000 (which they believed met the requirements of fully 
developed turbulent flow) . The profile was extrapolated to 
a universal profile at higher Reynolds numbers in an 
arbitrary manner due to lack of additional useable data. 

The chosen profile was required to satisfy continuity aud 
the Blasius friction factor. In the laminar sublayer the 
eddy diffusivity was taken as zero. In the turbulent core 
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a modified form of the mixing length expression presented 
Schlichting (13) was employed, with the coefficients 
determined ns functions of the hulk Reynolds number . The 
resulting expression for the velocity was quite complex and 
required numerical solution. Their predictions of fully 
developed Nusselt numbers for the low Reynolds number 
turbulent range showed a substantial improvement over the 
existing analyses since the latter had been based on the 
asymptotic velocity profile. 

lioelter, Martinelli, and Jonassen (l4j have also 
indicated that, for low Reynolds number turbulent flow, the 
distance to the edge of the laminar sublayer, y*, is not 
constant* This conclusion is further supported by the 
trend of the data obtained by Rothfus and Prengle (15) with 
a dye tracing technique to show the breakup of laminar flow. 
For the flow regime presently under considerat ion , the 
foregoing information invalidates the approximation that 
the laminar sublayer extends to a constant value of y + 
(commonly taken as 5) as utilized by some investigators 
(l(t, 16 ) for high Reynolds number flow. 

Adiabatic tube flow studies by Rotta (17) have 
shown int ermitt ent laminar and turbulent flow can exist in 
the transition flow regime. As the axial distance 
increases, the percentage of the time during which smooth 
flow is observed diminishes. Flow of this transitory 
nature, indicated by an intermit tency factor (the fraction 
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of time during which the flow .it a given position remains 
turbulent) of less than one, can persist for large axial 
distances before the development of a continuous fully 
developed turbulent flow (defined by Y ~ l). Or the 
transition may be completed in a very short distance. The 
apparent distance required shows a strong variation with 
Reynolds number in the range 2000 5= Re 3000 . Extrapola- 
tion of Rotta's results by McEligot, Ormand, and Perkins 
(9) yielded an intermi ttency factor of 0.95 at about 
twenty-five diameters for Re = 3000 , indicating fully 
turbulent flow can be achieved at this Reynolds number for 
tubes with reasonably short hydrodynamic entry lengths. On 
the other hand, Preston (l8) examined Nikuradse's 1932 and 
1933 friction data and noted that data for all roughnesses 
tended to a smooth tube friction curve, from which a 
transition to the laminar curve began to occur at a bulk 
Reynolds number of 4470. He concluded that this was the 
lowest possible Reynolds number for fully developed flow. 
However, Knudson and Katz (l), Senecal (4), and Schlichting 
(13) are a few of numerous authors who present data that 
are in agreement with the Blcisius smooth tube friction 
curve for Reynolds numbers as low as $000 , The apparent 
disagreement is probably explained by the observation that 
Nikuradse's friction data for Re « 5000 involved total 
pressure differences with magnitudes of order of the 
experimental uncertainty (19)* To avoid such conceptual 
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difficulties in the present paper, fully developed flow 
will be defined in the terms of the intermittency factor, 
i . e . , - 1 . 

Sources of velocity traverse data in the Reynolds 
number range of interest, 3000 to 15 * 000 , are presented in 
Table 1. It is seen that data are scarce. Except for the 
data of Senecal (k) and Nikuradse (19)» data are untabu- 
lated and are thus not available in a form of sufficient 
accuracy for development of a velocity profile. 

Nikuradse' s data (19) were collected in the exit plane of a 
circular tube where excessive free turbulence probably was 
being generated and transverse mean flow was not inhibited 
by the tube wall. It is concluded that only the work of 
Senecal--for the narrow range, 3000 Re <= 4000--are useful 
for the present purposes. Additional data are necessary to 
describe the approach from Re = kOOO to the asymptotic 
profile . 
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Table 1. Available Adiabatic Low Turbulent Reynolds Number 
Circular Tube Velocity Data for Re c 15 *000 


Invest i gators 

Re 

Data 

tabulated 

Tube ID, 
inches 

Loc a t ion of 
measuring 
station , 
diameters 

Nikuradse (19) 

4,000 

6,100 

9,200 

Y es 

0.3937 

0.001 to 
0.002 

behind exit 
pi aii e of 
tube , i . e . , 
outside 
tube 

Rothf us , Monrad, 

6,210 

No 

3.0 

« 100 

and Senecal (20) 

Seban and 

6 , 200 

No 

1.49 

90 

Shimazaki (21) 

Deissler (22) 

8,000 

11,000 

i4,ooo 

No 

r^ 

CO 

• 

0 

100 

Senecal (4)* 

3,002 

3,464 

4,108 

3,062 

4,085 

Yes 

0.50 

0.75 

360 

240 

Bak ewe 11 ( 23 ) 

8 , 700 

No 

11.2 

~ 27 


•Senecal presents additional profiles at Re «c 3000, 
but for these data it is not clear that ^ would be approxi- 
mately unity at the measuring station. 




CHAPTER III 


VELOCITY DISTRIBUTION 

A continuous, Reynolds number-dependent , veloci ty 
profile is developed. It is based on a modification of 
Reichardt’s wall and middle law eddy diffusivity expres- 
sions (2*t) and additional adiabatic velocity measurements 
for fully developed turbulent flow. 

Experimental Apparatus 

The facility employed for the measurement of the 
velocity and wall friction data was located at Texas 
Technological College. It has been described in detail 
elsewhere (25)*^ Important details of the apparatus, and 
the modifications performed for the present study are 
included herein. The basic test section was horizontal and 
unheated. It was fabricated from a twenty foot section of 
1.6l inch i.d. Inconel pipe. Twenty-seven static pressure 
taps (with l/l6 inch diameter holes) were axially spaced at 
2.720 inch intervals upstream and downstream from a probe 
traversing mechanism. The traversing mechanism was relo- 
cated to 255 diameters from the inlet by adding another 

1 . Copies available upon request from the Depart- 
ment of Mechanical Engineering, Texas Technological 
College, Lubbock, Texas. 
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twenty foot section ot' Inconel, pipe, with identical dimen- 
sions, upstream from the basic test section. Mating ends 
of the two pieces of tubing wore machined square and stool 
collars were used to assure alignment at the joints. 

Joints were made air tight by applying U. S. Royal indus- 
trial adhesive (epoxy) at the ends of the collars. 

Alignment of the two sections was achieved with a six foot 
carpenter's floor level. An exit length of 62 diameters 
followed the probe. Details are presented in Appendix A, 
Figure A-l. 

All velocity measurements were obtained with a 
square ended impact probe constructed from 0.0355 inch o.d. 
stainless steel hypodermic tubing. The particular probe 
used had an impact tube coefficient of one since calibra- 
tion showed that viscous corrections were unnecessary for 
local probe Reynolds numbers greater than 20 (see Appendix 
A). The particular probe design, used is also the least 
sensitive to pitch and yaw (25)* 

The probe traversing mechanism consisted of a 
micrometer drive, graduated in 0.001 inch increments, with 
an Inconel spindle. A O.O 36 inch hole had been drilled 
along the axis of the spindle and the velocity probe was 
silver soldered to the end. The probe traversing mechanism 
was housed within a 1/2 inch diameter thin wall tube, which 
was welded to the test section at a right angle to the wall 
static taps. A similar arrangement, use to obtain 
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temperature profiles, was located diametrically opposite 
the velocity traversing mechanism. A small horizontal slit 
was milled through the test section wall lor each mechanism 
to allow the probe to he completely withdrawn into its 
housing when not in use (see Appendix A, Figure A-2). In 
this construction the velocity probe tip was constrained to 
travel along the radius located at 90 ° from the probe 
reference wall static tap. The temperature probe was not 
utilized during the present investigation . 

A Flow Corporation Model MM3 Micromanometer, with a 
resolution of +0.0001 inches and a stated accuracy of 
+0.0002 inches of butyl alcohol, was employed to measure 
the impact (velocity) pressure drop and the axial wall 
pressure drops. The velocity was calculated from the 
expression 


-W 2g c AP 

P 


(3-1) 


The practice of neglecting the contribution of the 
fluctuating turbulent velocity components, which are of 
importance only in the immediate vicinity of the wall, is 
incorporated into the above equation ( 26 , 27) • Extrapola- 
tion of Sandborn's ( 28 ) fluctuating velocity component data 
to the low Reynolds number range, suggests the error 
created by this omission is probably negligible beyond one 
probe diameter (y/r = 0.044l) from the wall. For the 
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present investigation this distance corresponds to a y* 

4.3 to 12 corresponding to the lowest and highest Reynolds 
number runs, respectively. En his work, y/r = 0.01 was 
the largest value at which a probe correction was utilized. 
Hanratt> and Johnk (11 ) and Daily and Hardison (2?) indi- 
cate that disagreement and uncertainty exist in the 
literature concerning probe corrections. Although the work 
of Daily and Hardison contained an extensive and compre- 
hensive literature survey on impact probes and corrections, 
neither group used a correction for their data. The main 
justification in both cases was the inability to formulate 
meaningful, accurate correlations from the available 
literature . 

In the present work, the corrections suggested by 
Daily and Hardison, including Page's correction for the 
effective center displacement of a probe in a circular tube 
( 29 ) > have been considered with the uncertainty analysis. 
This analysis revealed the uncertainty in measuring the 
probe pressure drop always predominated over the possible 
error in neglecting the fluctuating velocity correction. 

At y/r values of 0.0472 and 0.1093 and a Reynolds number 
of 4o8o, the difference in the predicted uncertainties with 
and without consideration of the corrections are 3*64 
versus 3*43 and 2.64 versus 2.58 per cent respectively. 

The corrections and uncertainty analysis are discussed in 
Appendix A. 
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A 0-125 psig compressor supplied air througli a 
large storage tank. Three pressure regulators connected 
i»i series controlled the flow. Ihe last regulator was 
capable of controlling with 0.01 psig sensitivity and 
0.013 psi drift in 15 hours lor upstream pressures of 0 to 
23 psig. For flow at Reynolds numbers of 7000 and above, 
the last two regulators were bypassed. During the collec- 
tion of the friction data, the probe was completely with- 
drawn into its housing. 

Experimental Results 

Adiabatic velocity profile and friction data were 
collected for Reynolds numbers of 3020, 4o8o, 5010, 7030, 
10,100, and 15*000. Fully turbulent flow, as defined via 
the intermit tency factor, was r-hecked with hot wire 
measurements of the fluctuating axial velocity component 
at the probe station (255 diameters) and at approximately 
307 diameters (7*5 diameters from exit). The inter- 
mittently smooth and jagged output, observed at lower 
Reynolds numbers, was completely undetectable above a 
Reynolds number of 2700. 

Wall friction data are presented in Figure 1 for 
the 1.6 inch tube used in the velocity study and the 1/4 
inch tube used in the heat transfer study. For the 1.6 
inch tube, good agreement is seen between the Blasius 
friction factor expression (1) and the experimental data, 
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except tor the lowest KojuoMs number* run wheru the 
uncertainty is largest. 

The velocity data arc presented in standard non- 
dimensional coordinates (u + , y + ) in Figures 2 and 3* 

Shown, tor comparative purposes. are Scuecal ' s data for a 
bulk Reynolds number of 3002, as well as the semi- 
theoretical profiles to be discussed later. Uncertainty 
intervals for selectod data points, predicted after the 
manner of Kline and McClin took (30), are also presented 
In forming u + and y + the experimental friction data were 
employed. Dew point measurements showed that dry air 
properties were applicable. The test section temperature 
and pressure were considered invariant across any given 
cross section. A maximum Mach number of 0.023 was observed 
for the 15)000 Reynolds number run. 

As a check on the experimental velocity profiles, 
the bulk velocity was obtained by graphical integration 
(with a planimeter). The maximum percentage difference 
between the integrated and measured mass flow rates was 
+3*6 per cent at a bulk Reynolds number of 3017* For this 
Reynolds number, Fage (29) indicates an error of approxi- 
mately +1.4 per cent can exist in the integrated bulk 
velocity due to the effective center displacement effect. 

2. Details of the data collection and reduction 
are presented in Appendix A. Tabulated values are available 
in Appendix B. 
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An additional error is present due to the flow area reduc- 
ti on that occurs when the probe is inserted into the tube* 
With the probe fully extended a maximum reduction of 1*4 
per cent occurs in the flow area. The width of the milled 
slot in the wall is approximately 1 per cent of the 
circumf erence . 

A large deviation from accepted viscous sublayer 
magnitudes is observed in the immediate vicinity of the 
wall (y+ «e 15 ). The deviation can be attributed to the 
presence of the small milled slot in the tube wall. This 
disturbance in the wall could cause the generation of 
additional turbulence in the proximity of the slot resulting 
in values of the impact pressure drop differing from the 
pressure drop which would exist if the slot was not 
present. This deviation eliminated the possibility of 
using these data to determine effective laminar sublayer 
thicknesses and laminar sublayer velocities accurately. 

From a qualitative examination of the data, it appears the 
slot effect is no longer felt beyond y/r ranging from 
approximately 0.09 to 0.05 (y + * 9«7 to 21. 3) t correspond- 
ing to the lowest and highest flow rates, respectively. 
Experimental velocity profiles, presented by Haugen and 
Dhariak (31) for plane flow over rectangular slots, show 
that such slots may be expected to have a negligible 
effect at y + * 40. Unfortunately , their results were 
presented graphically and could not be used to establish 
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tlie lower bound for whic„. the present data would be valid. 
For further details the reader is referred to Appendix A. 

Semi -Empirical Profile 

Theoretical objections to the earlier McEligot, 
Ormand , and Perkins (9) velocity formulation exist, even 
though substantial improvement in the low Reynolds number 
heat transfer predictions was obtained with its use. Their 
development produces zero Vctlues for the eddy diffusivity 
in the laminar sublayer and at the tube centerline. 

Current thinking supports the concept of a smooth eddy 
diffusivity distribution which varies as the fourth power 
of y on approaching the wall (2, 3)» Moreover, their 
profile is based on the extrapolation of data in the 
narrow Reynolds number range between 3000 and 4000 . The 
present development was undertaken to resolve these objec- 
tions . 

Using his data as a basis, Reichardt (2*0 has shown 
that for high Reynolds number, fully turbulent, tube flow, 
the eddy diffusivity in the turbulent core can be repre- 
sented by a middle low which has a non zero value at the 
centerline 

V ■ ^ [‘ - * “t’l 


Near the wall he suggests the form 
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^ - #(y + - y* tnuh yVy*) (3-3) 

as reasonable, with the values H ~ 0.4 and y* = 11. In the 
present study a continuous expression was developed by 
expressing equation 3-2 in the equivalent form, 


- *Y* 

V ~ * 


2 - 


ll 


1 + 2(1 - 


+ 2 

*r) 


( 3 - 4 ) 


and by replacing #y* by the right side of equation 3-3 » 
which closely approaches Xy + in the core. For the result- 
ing expression, 



y^ tanh V 
y i 




1 + 2(1 



( 3 - 5 ) 


ft and y* were determined as functions of bulk Reynolds 
number, as described below. 

For fully developed flow in a circular tube, a force 
balance shows the shear stress distribution to vary 
linearly with the radius. Introduction of the defining 
equation for shear stress (equation 2-1) yields 


(1 - 


(1 




( 3 - 6 ) 


where yt = 7 p- ~\J t' (Re) /2 






2k 

The velocity is obtained directly by integration of this 
relation, with ^ m/\J described by equation 3-3* Numerical 
rosults cannot be obtained without assigning values to £ 
and y* . The bulk velocity may be defined as the one- 
dimensional velocity equivalent to the mass flow rate, 
i.e., for incompressible flow 



(3-7) 


The Blasius friction factor 

f = f (Re) » 0.0791 Re" 1 / 4 


(3-8) 


was chosen to evaluate the friction factor for the present 
study . 


From the definition of the friction factor it can 
be shown that the dimensionless bulk velocity and the 
friction factor are related by 



(3-9) 


For a given friction factor (Reynolds number), only a 
discrete number of combinations of # and y* exist which 
will satisfy continuity, i.e., produce the same value for 
the bulk velocity from both equations 3-7 and 3-9 • At 
various Reynolds numbers, a trial and error solution was 
performed on the University of Arizona IBM 7072 digital 
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computer to obtain possible combinations which would 
satisfy continuity. Simpson's Rule was employed in all 
integrations performed. 

For each experimental Reynolds number, the most 
acceptable combination of y* and X , along with a limiting 
band of combinations, was determined; this band was plotted 
against bulk Reynolds number. The criterion employed for 
an acceptable combination was close agreement between the 
data and the velocity profiles predicted by numerically 
integrating equation 3-6 (with the given combination of A 
and y*) . As further guide in the selection, eddy diffusi- 
vities predicted by equation 3-5 For the same combinations 
were compared with eddy dif fusivities derived from the 
experimental data via equation 3-6. The non-dimensional 
velocity derivative was evaluated by graphically determining 
the slopes of several plots of a given profile and 
averaging the results. The plot of the acceptable combina- 
tions revealed that a constant value of % s 0.4225 would 
suffice. Values of y* corresponding to the value of 
chosen were then plotted on logarithmic coordinates, and an 
equation was developed to fit the data. 

It was found that the results for y* and /f, as 
functions of bulk Reynolds number could be represented by 
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y* = 11 + 9.1116 oxp(-0.272'*9 Re x 1C” 3 ) 

+ [l5.83 exp( -O . 9^98 Re x lO" 3 )]^ 13-10) 

ft s 0.4225 

The second exponential terra in the expression for y* is 
negligible for Reynolds numbers above 4000 . The expres- 
sions for y* and At are valid for Reynolds numbers to 

5 

approximately 10 . The limiting upper Reynolds number may 
be determined by the user as the point at which he is no 
longer willing to accept the Blasius correlation for the 
friction factor. For high Reynolds number flow the 
profile converges to Reichardt’s form with y* * 11 and 
ft = 0.4225* The use of a value of ^ larger than used by 
Reichardt has been indicated by Hinze (32) , and more 
recently has been employed by Spalding. Spalding suggests 
that a value of ft as large as 0.43 is not unreasonable 
(33). 

The velocity predictions are compared with the 
experimental data for Reynolds numbers of 3020 , 4080 , and 
5010 in Figure 2, and 7030, 10,100, and 15,000 in Figure 3* 
For comparative purposes, the McEligot, Ormand, and 
Perkins profile and the universal logarithmic profile 
(1, 9) are also included. 

Due to the method of evaluating the Reynolds 
number-dependent coefficients of the present formulation, 


the value of the predicted friction factor, 


f 



(3-U) 


must agree closely with the Blasius value when integration 
of the theoretical velocity profile is performed to obtain 
V*. The same applies for the Mc-Eligot, Ormand , and 
Perkins profile (9) which was also forced to fit the 
Blasius friction factor expression. The extension of the 
asymptotic profiles predicts friction values differing by 
30 to 7 per cent from the Blasius values for the low 
Reynolds number turbulent flow range of 3000 to 13*000. 
These results have been presented for comparison with the 
experimental friction factors in Figure 1. 



CHAPTER IV 


DISCUSSION OF ADIABATIC VELOCITY PROFILE 

As illustrated in Figures 2 and 3> the present 
formulation is seen to correlate the valid experimental 
data to within 7*5 per cent. For Reynolds numbers below 
4o8o the formulation of McEligot, Ormand, and Perkins (9) 
predicts higher velocities in the so-called buff er layer 
than the present formulation and appears to give a 
slightly better correlation of the buffer layer data at 
Reynolds numbers of 3020 and kQ&O • The present treatment 
predicts higher values in the turbulent core than the 
former formulation; substanti<al improvement in the predic- 
tion is noted at 3020. Both formulations predict the 
overall trends of the data, and both satisfy continuity 
(i.e., realistic friction factor) due to the method of 
determining the Reynolds number dependent coefficients of 
the velocity profile. The present velocity data do not 
clearly discriminate between the two semi -empirical pre- 
dictions. Since both formulations were forced to satisfy 
the Blasius friction factor, the friction factor measure- 
ments cannot be used to discriminate. The question as to 
which best describes low Reynolds number turbulent flow 
remains to be tested by the heat transfer studies. 
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From a theoretical standpoint the present profile 
is more desirable. The eddy diffusivity is finite to the 
wall, while the McEligot, Ormand, and Perkins (9) develop- 
ment assumes a zero value in a laminar sublayer.. Moreover, 

in the vicinity of the wall the present analysis approaches 

3 

zero as y , which Reichardt has shown from continuity 
considerations to be the limiting behavior for isotropic 
turbulence (Eirod (jk) has shown the power should be » k 
for random turbulence). Both the eddy diffusivity and 
velocity profiles are continuous over the entire tube and 
are seen to have zero derivatives at the center line. The 
present analysis also predicts a positive finite, value of 
eddy diffusivity at the center line, which Reichardt (2*l) 
and Hinze ( 32 ) have verified as correct, rather than the 
zero or minus one value predicted by many of the other 
existing formulations. Although the center line behavior 
of the eddy diffusivity has little effect on the results of 
heat transfer analyses, where the velocity -eddy diffusivity 
formula lions have found their main application, it is 
desirable for the formulation to satisfy as many of the 
known conditions as possible. 

It is not altogether surprising that y* was found 
to be Reynolds number dependent, whereas a constant value 
of iC was sufficient to correlate the data in the low 
Reynolds number range. In his original paper, Reichardt 
indicated that y* was a measure of the thickness of the 
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viscous sub layers, which references (0» 1^, 13 » 32 ) have 
shown to vary with the Reynolds number. Conversely, 
characterizes an eddy diffusivity expression valid 
primarily in the turbulent core of the flow, where the 
effect of the wall decays rapidly with distance. It is 
reasonable to expect that the value of /C would be equally 
valid at low Reynolds numbers as long as the flow was fully 
turbulent in the core and not effected by the proximity of 
another surface. As discussed earlier, the use of a value 
of yf larger than Reichardt's value of i).k is not a complete 
innovation. Reichardt * in his original paper, noted that 
the accuracy of the available measurements were not 
sufficient for a close determination of 1C . The main dis- 
advantage of the present formulation is that xt involves a 
numerical integration for the velocity profile. This is 
true of most of the more recent eddy diffusivity and mixing 
length models for both high and low Reynolds number 
turbulent flow. A closed form solution for the velocity 
cannot be obtained without additional simplifying assump- 
tions and relaxation of boundary conditions above and 
beyond those employed in the more recent works# 


CHAPTER V 


PREVIOUS DIABATIC STUDIES 

As a background to tho low Reynolds number turbu- 
lent thermal entrance region study, the present state of 
the art will now be considered. Experimental and 
analytical investigations of local heat transfer charac- 
teristics for turbulent flow in circular tubes, specifically 
directed at the Reynolds number regime from 3»00l* to 30,000, 
are almost nonexistent. Modern advances in technology have 
made it increasingly important for the design engineer to 
accurately' predict these characteristics for optimization 
within material limitations. Once thermal entrance region 
solutions are obtained for the constant wall heat flux 
boundary condition, it is possible to predict wall tempera- 
tures with axially varying wall heat flux by superposition 
(2). However, existing analytical formulations are 
inadequate for such predictions since they are based on 
high Reynolds number data and assumptions. This is 
especially true below bulk Reynolds numbers of 15,000. 

Analyses for the prediction of turbulent heat 
transfer and friction characteristics for fully established 
conditions and constant fluid properties are well estab- 
lished for Reynolds numbers in excess of 50 , 000 . The 
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analyses of Martinelli (lf>), Sparrow* Hallman, and Siegel 
(33)« l>r.\ssler (3fi), and the Dittus-Hoelter experimental 
correlation (2), with a 0.02] coefTiciont, all correlate 
the fully developed data to within 6 per cent. In light of 
entrance data collected in the last ton years (37), the 
thermal entrance region solution of Sparrow, Hallman, and 
Siegel has gained wide acceptance over the boundary layer 
integral analysis performed by Deissler ( 36 ) and later 
extended by Wolf (3^)* The Sparrow, Hallman, find Siegel, 
solution is based on a slightly modified version of 
Deissler 's velocity and eddy diffusivity expressions (12) 
and the assumed equivalence between the eddy diff usivities 
of heat and momentum. The form of the partial differential 
energy equation used was linear because fluid properties 
were taken to be constant. Since their treatment forms the 
basis for the present analysis it is described now in 
greater detail. 

Sparrow and co-workers separate ! the temperature 
field into two parts which represent the entrance and 
fully f developed regimes, respectively. An ordinary dif- 
ferential equation was obtained with the temperature in 
the fully established region as the dependent variable. 

This equation was integrated directly. For the remaining 
portion of the temperature field, a second separation of 
variables yielded (l) an ordinary differential equation 
which has an exponential solution for the axial variation, 
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and (2) a Sturin-Liouville equation for tho radial depend- 
ency. For the solution of the Sturin-Liouville equation, 
tho first five eigenvalues wore obtained at Reynolds 
numbers of 5^*000* 100,000, and J00,000 by employing a 
digital computer. Tho numerical mothod was not indicated, 
but the nature of the problem suggests a trial and error 
solution was probably used to obtain the. eigenvalues. 

In addition to the large number of constant 
property investigations , several recent variable property 
thermal entrance and downstream studies have been conducted 
(39* 40 , 4l ) . These studies have net been completely 
successful in predicting variable property uata. 

In the portion of the low Reynolds number range 
from 3,000 to 15*000 heat transfer data for gases are 
scarce and erratic. In the remaining range from 15*000 to 
30,000 data are more plentiful, but still not abundant . 

The majority of the data are confined to fully developed 
conditions or to average values for the entire length of a 
test section. McEligot, Ormand, and Perkins (9) point out 
an increasing deviation between the fully developed data 
and previous analytical predictions with decreasing 
Reynolds number. The discrepancy approached fifty per 
cent. This trend was also present in the majority of tho 
variable property data they examined. Molina and Sparrow 
(42) have presented variable and constant property entrance 
region data for water and oil (Pr = 3* 48, 75)* with the 
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Reynolds number range betweort 12,000 and 30,000 included. 
Their entrance region data show the dimensionl ess entrance 
distance (defined as the axial distance for the local heat 
transfer coefficient to reach ninety-five per cent of its 
fully developed value) to decrease with increasing Reynolds 
numbers. The same trend is indicated by the data of 
Hanratty and Johnk (11). In contrast, the analysis of 
Sparrow, Hallman, and Siegel (35) slows a very slight 
increase in the entrance length for increasing Reynolds 
numbers. The lowest Reynolds number considered in their 
solution was 50 , 000 . 

The increase in velocity with respect to the value 
predicted by a universal velocity profile (see Figure 2) 
can explain the deviation between heat transfer analyses 
and data for Reynolds numbers below 30,000. The observed 
thickening of the viscous sublayer represents a reduction 
in the turbulent component for momentum transfer. There- 
fore, the turbulent transport of energy is reduced so that 
the resulting Nusseit number decreases also. The available 
analyses, with the exceptions to be discussed below, are 

p 

formulated using semi-empirical velocity and eddy dif- 
fusivity, or mixing length, expressions based on universal 
velocity profiles which do not vary with Reynolds number. 

It should be expected that these analyses would be valid 
only for high Reynolds numbers. It may be recalled that in 
the wall frxution problem the comparable analyses showed 
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agreement with friction factor measurements in the high 
Reynold* number range only. 

Heat transfer analyses which treat the deviation 
from the universal velocity profile are scarce. McEligot, 
Ormand , and Perkins (9) predict Nusselt numbers, for fully 
established temperature fields in circular tubes, which are 
in agreement with available data. Their analysis utilized 
a varying velocity profile which was based on Senecal's 
data (4) at Reynolds numbers of 3»000 and * 4,000 and was 
arbitrarily extrapolated to a universal velocity profile at 
higher Reynolds numbers. Equivalence between the eddy 
dif fusivities of heat and momentum was employed. The 
energy equation was solved by numerical means. The Dittus- 
Boelter correlation with a coefficient of 0.021 correlated 
their data, and the data of others, to within five per cent 
of their predictions for Reynolds numbers above 4 , 000 . 

Gill and Lee (43) used the low Reynolds number 

velocity profile of Gill and Scher (7) to predict neat 

* 

transfer results for parallel plates with constant wall 
temperatures . McEligot, Ormand, and Perkins (9) point out 
that the results show a sharper approach to the laminar 
Nussext number than is noted in circular tube experiments. 
Haber&troh and Baldwin (44) adopted Pai ' s low Reynolds 
number profile (10) to predict temperature profiles and 
Nusselt numbers for fully developed tube flow with constant 
wall heat flux. For low turbulent Reynolds numbers, the 
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validity of this analysis, as well as the formulation of 
Gill and Lee (^3)> is questionahl e . Both studies were 
based on velocity profiles with monotonir. variation from 
high Reynolds number results to parabolic laminar profiles. 
Such a procedure is in disagreement with Rotta’s (17) 
experimental observations of intermittent laminar and 
turbulent flow at values of Reynolds numbers in the transi- 
tion range. 

Temperature profile measurements by Hanratty and 
Johnk (11) indicate the eddy diffusivity of heat varies 
with axial position in the thermal entry. These changes 
are confined to the centerline region of the tube after the 
first three to four diameters. But, they conclude the 
calculation method of Sparrow, Hallman, and Siegel (35) » 
coupled with Reynolds analogy, = £* , is adequate for the 
prediction of heat transfer results in the thermal entry 
region. Deissler (k$) Bad shown earlier that extension of 
adiabatic profiles was valid for moderate heating rates. 

The lowest Reynolds number employed in the work of Hanratty 
and Johnk was 18,000, but only the data collected for 
Reynolds numbers of 2^1,000 and 35*000 were sufficiently 
accurate to analyze the variation of eddy diffusivity in 
the entrance region. 

Of the limited nuitobuf of low Reynolds number 
adiabatic velocity profiles, only the formulation dua to 
McEligot, Ormand, and Perkins (9)» and the present profile 
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developed in Chapter 111 appear to be without serious 
objections. Though the present formulation is theoreti- 
cally more rigorous, both apparently correlate the low 
Reynolds number velocity data reasonably well. A substan- 
tial improvement in the prediction of fully developed 
Nusselt numbers has been obtained by McEligot, Ormarid , and 
Perkins (9) employing their profile. 

As noted, no thermal entry solution is available 
for low Reynolds number turbulent flow. The method of 
Sparrow, Hallman, and Siegel (33) will be applied in the 
following study, with both of the above formulations, (a) 
to provide a usable low Reynolds number thermal entrance 
solution and (b) to resolve, if possible, which formulation 
provides the best description of the velocity profile. 



CHAPTER VI 


HEAT TRANSFER ANALYSIS 


For steady, low Mach number, hydrodynamically fully 
developed flow of an incompressible fluid, with constant 
properties, negligible viscous dissipation and negligible 
natural convection effects, the energy equation for 
circular tubes can be written as 
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For the nondimensional form of equation 6-1, 
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a solution for the temperature profile and Nusselt number 
was obtained tit Pr = 0.71 for all axial positions. The 
boundary conditions were constant wall heat flux and 
uniform initial temperature. The method used by Sparrow, 
Hallman, and Siegel (35) was employed to solve equation 6-2 
with the following important changes: 
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1. Tho eddy diffusivity and velocity profiles for 
adiabatic low Reynolds number flow were used in 
place of tho profilo used by Sparrow, Hallman, and 
Siegel. The use of Reynolds analogy, - £ , was 
retained. It leads to 



with ft - 0.4225 



tanh 





(6-3) 


and y* = 11 + 9*1116 exp (-0.27249 K© r 10~ 3 ) 
+ [l5.83 exp (-0.9498 Re x 10~ 3 )1 


for the eddy diffusivity of heat. The veloci ty is 


obtained by integration of the defining equation 


for 


V ? 



(6-4) 


These profiles were shown to be valid for all fully 
turbulent Reynolds numbers greater than, or equal 
to 3000 , whereas the profiles used by Sparrow, 
Hallman, and Siegel are valid only above 30 ,000. 

2. The distance from the tube wall, y, was used in 

place of the radial coordinate as the independent 


variable 


3. A different nondimensiotializiition was employed. 

This leads to different numerical magnitudes for 
the eigenvalues than would be obtained using the 
nondimensionalization of reference (35) • 

Results were obtained numerically on an IBM 707^5 a 
fifth order Predictor-Corrector, with a Runge-Kutta proce- 
dure for starting, was employed in the eigenvalue solution. 
Since the series solution was truncated, results are not 
available at x/D = 0. Details of the numerical work are 
available in Appendix G. 

For later comparison, the same numerical program 
was used to solve the energy equation with the eddy dif- 
fusivity and velocity profile expressions used by Sparrow, 
Hallman, and Siegel (35) » and the low Reynolds number 
expression developed by McEligot, Ormand, and Perkins (9). 

A solution employing the latter profile was deemed desirabl 
since the velocity data presented in Chapter III did not 
discriminate between their profile and the present one. 

Both low Reynolds number profiles predict the deviation of 
the nondimensionalized veloc ity, u+, from the universal 
velocity profiles for Reynolds numbers below 30*000. This 
is not true of the profile used by Sparrow, Hallman, and 
Siegel . 

Due to the wide spread knowledge and availability 
of the method of Sparrow, Hallman, and Siegel (35) (see 
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Kays (2) and summary in Chapter V) , only the results 
obtained using equations 6-3 and 6-4 for equation 6- 2 are 
presented. Detuils of the present solution may be found in 
Appendix H. Comparisons with the solutions obtained using 
the other two profiles and data are deferred to a later 
section . 


The local A'usselt number may be given in terms of 
the fully developed value as 


Xu 

‘ Nl bo 



(6-5) 


The first seven values of the eigenvalues,^^ and A^, are 
presented for gases in Tables 2 and 3 for a range of 
Reynolds numbers. Since a searching method was employed in 
their evaluation, the accuracy of any individual eigenvalue 
is independent of previously calculated values. Eigen- 
values and associated coefficients for laminar flow, as 
calculated by the present numerical program, are also 
compared with the Siegel, Sparrow, and Hallman (46) laminar 
solution to establish confidence in the validity of the 
present numerical method. For the laminar case, :o 
ambiguity in the calculation of the velocity profile is 
possible. At the upper end of the range, Re = 50 , 000 , the 
Sparrow, Hallman, and Siegel solution (as generated by the 


Table 2. Eigenvalues 



o 0.02053 0.06705 0.13926 0.23709 0.36056 0.50963 


able 3* Values of 
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present program) is a maximum of 1.1 per cent higher than 
the solution obtained from their published eigenvalues. 

The Nusselt numbers for fully established condi- 
tions , used to normalize the thermal entry values, are 
represented within 4 per cent by a Dittus-Boel ter equation 
of the form 

Nu = 0.021 Re* )# ^ Pr^ * * ( 6 - 6 ) 

00 

for the range of Reynolds numbers from 3000 to 30 , 000 . The 
agreement could be improved to 3 por cent with the use of a 
coefficient of 0.0231 and a power of 0.79 for Re, but such 
reduction is probablv not worthwhile for design. 

The entry region solution is presented graphically' 
in Figure 4 and is listed in Tabie 4. The analysis shows 
the dimensionless entrance distance to decrease slightly 
with increasing Reynolds number. For all Reynolds numbers 
the Nusselt number can conservatively be considered to 
reach its downstream value by x/D = 30 . In the entrance 
region the Nusselt number increases markedly, at a fixed 
axial position, as the Reynolds number decreases. The main 
effect is confined to distances less than 10 diameters; at 
x/D * 3 , a l 6 per cent increase in the normalized Nusselt 
number is observed between Reynolds numbers of 50,000 and 
3,000. The solution is approaching the laminar solution of 
Siegel, Sparrow, and Hallman (46), though a considerable 










Table 4. Entrance Region Nusselt Number 
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Re 

3,000 

4,000 

6,000 

10,000 

20,000 

50,000 

Nu oo 

11.056 

l4 .081 

19.571 

29.463 

50.275 

100.810 

v /n 
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2 • 5 

1.4951 

1.434 

1.3780 

1.3292 

1.2936 

I.2638 

5.0 

1.2655 

1.2306 

1 . 2020 

1.1795 

1.1657 

1 .1550 

7.3 

1.1674 

1.1445 

1.1273 

1.1148 

1.1087 

1.1049 

10.0 

1.1126 

1.0967 

I.O856 

1 .0782 

1.0757 

1.0751 

12.5 

1.0784 

1.0669 

1.0595 

1.0549 

1.0543 

1.0553 

15.0 

1.0536 

1.0471 

1 .0421 

1 .0392 

1.0397 

1 .o4l5 

17.5 

1.0399 

1.0335 

1.0301 

1 .0783 

1.0292 

1 .0314 

20.0 

1.0288 

1 .0240 

1.0216 

1 .0205 

1 .0217 

1 .0240 

22.5 

1 .0209 

1.0173 

1.0156 

1 .0150 

1 .0162 

l .0184 

25.0 

1.0152 

1 .0125 

1 .0113 

1 .0109 

1.0121 

1 .0141 

27.5 

1 .0111 

1.0090 

1 .0082 

1 .0080 

1 .0090 

l .0109 

30.0 

1.0081 

1 .0065 

1.0059 

1 .0059 

1 .0068 

l .0084 

32.5 

1.0059 

1 .0047 

1 .0043 

1.0043 

1.0051 

l .0065 

35.0 

1.0043 

1 .0034 

1 .0031 

1.0031 

1.0038 

1.0050 

37.5 

1.0032 

1.0025 

1.0023 

1.0023 

1.0029 

1.0039 

4o.o 

1.0023 




1.0022 

1.0030 
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difference still exists in the normalized Nusselt number 
between the turbulent values at 3000 and the laminar 
results at 2000. Doth the Nusselt number and th~ entry 
length variations are consistent with the trends of the 
laminar solution, the oil and water data of Malina am, 
Sparrow (42), and the air data of Hanratty and Johnk (ll). 
At a Reynolds number of 30,000 the solution, which is based 
on the present velocity profile, agrees within two per cent 
with the solution of Sparrow, Hallman, and Siegel as 
calculated by the present program. Agreement to within one 
per cent is obtained using their published eigenvalues. 

To ease application by the engineer, the entrance 
region results for the constant heat flux boundary have 
been correlated by the equation 

- 2 -1 

§* 7 - * 1 + 0.8(1 + 70,000 Re 2 )(j£) (6-7) 

wu oo ' 

Agreement is within 5 per cent for x/D a* 2, with improve- 
ment to 3 per cent or less for x/D =* 4 over the range 
3000 < Re < 50,000. For x/D s» 12 the prediction could be 

JjJ n 

improved by the use of a separate correlation for 
similar to the above expression. The somewhat "poor” 
disagreement (3 per cent) in the downstream region is 
tolerated for the sake of simplicity in the correlation. 

The analysis shows that the fully developed value is 
reached for x/D 30. 
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To treat problems with axially varying wall heat 
flux the local Nusselt numbers are not directly available. 
Instead, they must be derived for the heat flux variation 
of interest by employing the eigenvalues and constants. 

The eigenvalues and associated constants of the present 
study were found to be correlated by the expi cessions 


A 


2 

n 


-b 


A, Re 
1 » n 


1 ,n 


-d 


+ C, Re 
1 ,n 


l,n 


( 6 - 8 ) 


and 

An = -A„ Re 2,11 + C 0 Re 2,11 (6-9) 

2 , n 2 , n 

for Reynolds numbers from 4000 to *K>»000. The above three 
correlations are equally valid for helium (Pr » 0.664) and 
air (Pr » 0.71). Numerical values for the coefficients and 
constants are pres ited in Table 5. 
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Table 5* Numerical Coefficients and Powers for Eigenvalue 
Correlation 


A_ 

l,n 

0.170 

b 

1 

0.1098 

C 

1.347 

1 ,n 
x 10^ 

d, 

1 , n 

2.5 

0.491 

0.1064 

1.466 

x 10*? 

2.34 

0 . 945 

0.1029 

2.39 

x 10^ 

2.33 

1.531 

0.1001 

2.65 

x 105 

2.30 

2.29 

0.1001 

3.32 

x 10? 

2.64 

0*33 

0.1033 

O.059 

x 10 5 

1.662 

^2,n 

1.25 

b 

2 ,n 
0.179 

c 

-2.84 

2 ,n 
x 10 

d 2,n 

3.64 

0.650 

0.172 

-5.91 

x 10 * 0 
x 10, 

X io 4 

1.061 

0.459 

0.171 

2.82 

3.99 

0.235 

0.1299 

9.86 

1.885 

0.1193 

0.0817 

1.25 

x 105 

1.91 

0.0442 

0 

5.13 

x lo5 

2.13 
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CHAPTER VII 

EXPERIMENTAL HEAT TRANSFER INVESTIGATION 

Apparatus 

The apparatus used was a redesigned version of the 
h#at transfer loops used by Perkins and Worsoe-Schmidt (40) 
and McEligot, Magee, and Leppert (4l). Briefly, the test 
section consisted of a vertical, thin walled, one-quarter 
inch, circular Inconel 600 tube with a 100 diameter hydro- 
dynamic entrance, 100 diameter heated length, and a 30 
diameter exit. Two pressure taps were located just beyond 
the extremes of the 100 diameter heated length. The test 
section was enclosed in a six foot long, six inch diameter 
vacuum chamber to minimize heat loss effects and allow 
localized heat loss calibration. The vacuum environment 
also reduces the response time necessary to reach steady 
conditions with heating. A. C. power for electrical 
heating was provided by a line voltage stabilizer, a 
variable transformer and a 20-to-l transformer in series. 
Above approximately 80 amperes, a Lincoln model TM-500 
AC/DC welder was used as the power source. Variable area 
flow meters, with specified accuracy of one per cent of 
full scale, were used as the primary flow .n^rnr Client 
devices. A Foxboro integral orifice, differential pressure 
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cell served as a secondary check. The gas was supplied by 
a four bottle gas manifold. Test gas passed through two 
pressure regulators connected in series and through the 
flow measuring equipment before entering the test section. 
Flow control was achieved by means of two regulating 
valves, located downstream of the test section. 

Method of Testing 

The method of testing closely parallels that of the 
above references ( 40 , 4l ) . In summary, the results 
presented herein were obtained for the flow of helium and 
air through a vertical, resistively heated tube. Premium 
grade Chromel-Alumel thermocouples were spot welded to the 
tube at various axial positions. Heat transfer data were 
obtained for a number of different low power settings at 
each of four different inlet Reynolds numbers in the 
desired range. Measurements of the tube wall temperatures 
were obtained for axial positions ranging from x/D = 1.21 
to 95*68 with a maximum wal \ -to-bulk temperature of 1 . 4l . 
Test section pressures of approximately 100 psia were 
employed to hold the entering Mach number below 0.033 to 
avoid compressibility and viscous dissipation effects. The 
modified Grashof /Reynolds quotient was kept below 20 to 
avoid natural convection effects (47). The reader will 
recall that these effects were neglected in the analysis to 
which the data are compared. To determine the thermocouple 
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conduction error of the wall thermocouples ( 48 ), internal 
test section temperatures were measured with a ceramic 
thermocouple probe when the system was at power but without 
flow; comparison to the outside wall thermocouples led to 
calibration of the error for each wall thermocouple. The 
test section emissivity, e(x, T), was also determined 
locally from these calibration runs. This calibration pro- 
vided a convenient method for evaluating the local external 
heat loss. 

The local energy generation was evaluated from the 

2 

i R' product. In the resistance calibration the assumption 
was made that the local resistance at a given temperature 
was equal to the average test section resistance at the 
same temperature. The variation of the average test 
section resist nee with temperature was obtained from 
voltage anl current measurements taken during a large 
number of runs without flow. For heating without flow the 
axial temperature remained essentially constant. These 
runs were spaced over a two year period and indicated a 
slight decrease in resistance with time. Consequently, for 
the resistance variation with temperature, the data reduc- 
tion procedure used a straight line correlation which 
favored the resistance data obtained during the same time 
period as the heat transfer data. The wall heat flux was 
evaluated by subtracting the radial radiation and axial 
conduction heat losses from the local energy generation. 


33 


Data reduction was performed on an IBM 7072 digital 
computer employing a modification of the data reduction 
program used by McEligot (^9) . The reduced data are 
presented in Appendix II. 

A method similar to that of Molina and Sparrow (^2) 
was used to approach the constant property idealization of 
the analysis. For each of the fixed inlet Reynolds 
numbers , the data for the various power settings were 
normalized by a Di ttus-Uoelter correlation (with a 0.021 
coefficient) and were plotted against the wall-to-bulk 
temperature ratio on logarithmic coordinates . This choice 
of normalization partially removes the effect of small 
variations in the inlet Reynolds number, even though these 
variations never exceeded 0.8 per cent for a given flow 
setting. For a given thermocouple, extrapolation of the 
data to a wall-to-bulk temperature ratio of one removes the 
temperature dependent variable properties effect. The 
extrapolation procedure is demonstrated in Figure 5 for 
several of the thermocouples • 

For the helium data the uncertainty of the extra- 
polated constant property Nusselt numbers is estimated to 
range from a low of four per cent to a high of nine per 
cent for each of the inlet bulk Reynolds numbers of klBo f 
6800, and 10 , 300 . Typically, the uncertainty at first 
decreases with axial position and, then, begins to increase 
at large x/D due to increasing uncertainty in the bulk 
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temperature. A maximum difference of 0.8 per cent was 
observed between the mass flow rates measured witn the 
primary and secondary flow meters. Adiabatic friction 
factor measurements) taken before and after the heating 
runs without changing the flow settings) showed a maximum 
deviation from the Blasius friction factor (l) of -2 per 
cent; it occurred for the 10)300 Reynolds number runs. 

The uncertainty for the air data at a Reynolds 
number of 3^60 was considerably larger than for the helium 
Nusselt numbers; it ranged from 8-1/2 to 13 per cent. This 
increase is predominantly due to the very small range of 
wall-to-bulk temperature ratios for which the air data were 
obtained (a maximum of 1.09 versus 1 . 4 for the helium 
data). The percentage uncertainties in Nu are high for 
small wall-to-bulk temperature differences. Without data 
at larger ratios t the estimated uncertainty bands obtained 
in the constant property extrapolation process become 
larger. In addition) an increased uncertainty in the mass 
flow rate (5*7 per cent versus a maximum of 3*2 per cent 
for the helium runs) also contributed to the higher 
uncertainties for the air runs. A difference of 3-3/4 per 
cent between the primary and secondary flow rates and a 3 
per cent overprediction of the adiabatic friction factor 
were observed. This indicates the measured flow rate could 
be low, which would result in the calculation of an exces- 
sive value for the Nusselt number. 
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In summary, the dominant uncertainty in the Nusselt 
number was found to be due to the uncertainty in the wall - 
to-bulk temperature difference, which decreased in 
importance as the difference increased. However, most of 
the results can probably be considered conservatively to 
have uncertainties within 10 per cent. The reader should 
be careful not to be misled by the exaggerated expansion 
of the ordinate on some of the figures. 
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CHAPTER VIII 

COMPARISONS OF HEAT TRANSFER RESULTS 


Entrance Region 

The constant wail heat flux analysis is compared 
with the extrapolated constant property data for air at a 
bulk Reynolds number of 3 & 60 * and for helium at Reynolds 
numbers of 4l8o, 6800 , and 10,250 in Figures 6 through 9« 
The numerical solution was rerun for each set of experi- 
mental conditions to obtain the theoretical predictions 
shown. Solutions generated using the McEligot, Orraand , and 
Perkins (9)» and the Sparrow, Hallman, and Siegel (35) 
velocity profiles are also included. The idealization of 

constant wall heat flux is satisfied for x/D =* 2, with 
q” 

■y — having reached 0.955 to O .985 at x/D = "' l . 2 for all 
^ w max 

runs. The experimental heat flux distribution may be 
represented by an exponential variation in the immediate 
thermal entry, 


l" 

l w 


s q*» (l - 
4 wmax 


e 


-mx 


For this axial distribution the method of superposition 
(2) yields 
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Figure 6. Thermal Entrance Solutions Compared with Experiment 





Pr = 0.66 



Figure 7* Thermal Entrance Solutions Compared with Experiment 





Figure 8. Thermal Entrance Solutions Compared with Experiment 
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This solution approaches the constant wall boat flux 
solution to within two per cent or less by x/l> * k under 
the conditions of the present experiments. 

The air data are observed to be consistently higher 
than the analysis, a trend which is consistent with the 
suspected low value of the measured flow rate mentioned 
earlier. Excluding these data and the values at x/D = 

1 . 21 , which are highly inaccurate due to axial conduction 
in the test section wall and to the small wall-to-bulk 
temperature differences which exist near the initiation of 
heating, the data are correlated to within five per cent 
by the present analysis. Including all data, correlation 
is obtained to within ten per cent. 

From the present investigation , it cannot be 
definitely concluded whether a solution based on the 
present low Reynolds number velocity formulation, or on 
the McEligot, Ormand , and Perkins (9) formulation, is more 
accurate. Both correlate the data within the experimental 
uncertainty, and agree closely with each other, differing 
by a maximum of 6.7 per cent at a Reynolds number of 3000. 
Both satisfy continuity and converge to the Sparrow, 
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Hallman, and Siegel ( 33 ) solution at a Reynolds number of 
50*000 as shown in Figure 10. A slightly better agreement 
is noted for the present solution (2 versus 4-3/4 per cont 
maximum difference). The present formulation is recom- 
mended over the McEligot* Ormond, and Perkins solution (9) 
since it predicts lower* and thus more conservative, Nusselt 
numbers at the low end of the Reynolds number range. In 
addition* the eddy iiffusivity-velocity profile used is 
based on data covering the entire low Reynolds number 
range from 3000 to 15*000 and is theoretically more 
rigorous. The eddy diffusivity distribution satisfies the 
accepted requirements at the wall and centerline of the 
tube more completely than does the distribution of the 
earlier study. 


Fully Developed Solution 

The analytical prediction of Nusselt numbers for 
fully established conditions are compared with local data 
of the present work* and with the data presented by 
McEligot* Ormand* and Perkins (9)* in Figure 11. Both of 
the low Reynolds number solutions correlate the data to 
within a few per cent and represent a substantial improve- 
ment over the analyses based on universal velocity type 
distributions. As mentioned earlier, the fully developed 
Nusselt numbers of the present analysis are approximated 
within 4 per cent by the Dittus-Boelter equation 



Figure 10. Thermal Entrance Solutions at Re = 50,000 





Figure 11. Heat Transfer in Fully Established Flow: Semi -Empirical Predictions 
Compared with Eat a 
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Nu * 0.021 Re 0,8 Pr° * ^ (0-1) 

OO 

for the rang© of Reynolds numbers from 3000 to 50,000. 

The present formulation predicts a slight! lower 
Nusselt number than the McEligot, Ormond, and Perkins (9) 
solution at the lower Reynold*, numbers with the converse 
occurring above a Reynolds number of approximately 20,000. 




CHAPTER IX 


DISCUSSION OF HEAT TRANSFER STUDY 

The success of the present entrance region solution 
for the low Reynolds number turbulent heat transfer yic ds 
additional support to the use of Reynolds analogy, £* » 

in heat transfer analysis for gases. The data were 
correlated reasonably well by both of the low Reynolds 
number eddy 'if fusivity-velocity profiles considered, 
although they differed in several respects. The differences 
in the profiles indicate that a velocity profile is 
adequate for heat transfer predictions provided (a) it 
satisfies continuity, i.e«, predicted friction factors 
agree with data, and (b) it agrees approximately with the 
observed deviation of the velocity data from the universal 
profile • 

The eigenvalues of ths present analysis provide the 
necessary building blocks for the prediction of heat 
transfer results for axially varying wall heat flux at low 
heating rates. For the treatment of the method of super- 
position presented by Kays (2) the nomenclature are related 

a *7m 'An RePr “ nd A m ' ‘7m A n /N V whore An> A n> and Nu OO 

are the values predicted by the present solution. Care 
must be exercised if alternate references are used, since 
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the numerical values of the eigenvalues and the associated 
constants will differ as the nondimeusionalization differs. 

The present analysis is not necessarily restricted 
to Reynolds numbers below 50*000 since the semi -empirical 
velocity formulation used is valid for all Reynolds numbers 
for which the Blasius friction factor is acceptable* 

Results were not presented beyond 50,000, since the 
designer will probably wish to use the more familiar 
Sparrow, Hallman, and Siegel (35) solution at the higher 
Reynolds numbers. Below a Reynolds number of 50 , 000 , the 
present solution predicts slight increases in the entrance 
length and strong increases in the normalized Nusselt 
number in the entrance region with decreasing Reynolds 
numbers. The latter effect is most prominent for distances 
less than ten diameters. For a Reynolds number decrease 
from 50,000 to 3000 , increases of 3*5 per cent and l 6 per 
cent are observed in the normalized Nusselt number for 
x/D * 10 and x/D * 3* respectively. In all cases, the 
Nusselt number may be considered fully developed by 
x/D = 30 . As previously indicated, the above trends are 
consistent with the data of Malina and Sparrow (42), 
Hanratty and Johnk (ll), and the present data. Above 
Reynolds numbers of 50,000, a reversal in the trends of 
the solution was observed. Very slight increases in both 
the entrance length and the normalized Nusselt numbers 
occurred with increasing Reynolds numbers. This reversal 
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agrees with the small increases shown by the solution of 
Sparrow, llallman, and Siegel (35) for their published 
values of Reynolds numbers of 50,000, 100,000, and 500 , 000 . 
The change in the direction of the trends at high Reynolds 
numbers may be due to the reduced thickness of the laminar 
sublayers, and/or numerical inaccuracies in the solution. 

The reader's attention is also called to the 
requirement of fully turbulent flow, as defined by an 
intermit tency factor equal to one. The velocity profile 
employed is valid only for this condition. It is incapable 
of describing intermittent laminar and turbulent flow, 
which exists for values of the intermittency factor less 
than one. Thus, these results should not be used for 
prediction during the axial transition from laminar to 
fully-turbuleiit flow. 

For x/D values greater than four, the experimental 
Nusselt numbers were observed to decrease slightly as the 
wall-to-bulk temperature ratio increased from approximately 
1 to 1.4. These data might be correlated by an equation of 
the form 

^ c\ 

Nu « 0.021 Re 0,8 Pr° * 4 (=~) (9-1) 

A b 

but experiment'll uncertainties prevent an accurate evalua- 
tion of the exponent, a, since the range of T /T^ is small. 
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The question--which of the two low Reynolds number, 
turbulent eddy diffusivity and velocity profile formula- 
tions considered best describes low Reynolds number 
turbulent f low--remains unresolved. The differences 
between the two profiles noted in Chapter IV, however, 
appear to have a minor effect on the prediction of the heat 
transfer parameters, so resolution is probably not neces- 
sary . 

The use of the high Reynolds number expressions 
employed by Sparrow, Hallman, and Siegel (35) proved 
totally ineffective for thermal entrance region predictions 
in the low Reynolds number range. This is due to the use 
of a velocity formulation which neither satisfies 
continuity, nor describes the adiabatic profile accurately 
at low Reynolds numbers. For high Reynolds numbers, these 
requirements are fulfilled and the solution provides an 
accurate prediction of heat transfer characteristics for 
Reynolds numbers over 50,000. The agreement at high 
Reynolds numbers further confirms that adequate heat 
transfer predictions are possible only if both of these 
requirements are met for the velocity profile. 






CHAPTER X 
CONCLUSIONS 

The lack of accurate predictions of thermal entry 
heat transfer characteristics for low Reynolds number, 
constant property, turbulent flow has been eliminated. 
Design equations have been developed from this study to 
provide the engineer with a convenient method of analysis. 
The following conclusions can be drawn: 

1. A universal velocity profile does not provide an 
adequate prediction of the friction factor for low 
Reynolds number flow. It satisfies continuity only 
for Reynolds numbers in excess of approximately 
30,000. 

2. To be valid, a velocity formulation must predict 
both friction factor and velocity magnitudes 
consistent with the experimental datt\. The present 
modified form of Reichardt's eddy diffusivity 
yields a velocity profile which meets these 
requirements in both the low and high Reynolds 
number ranges. It is valid for all Reynolds 
numbers for which the Blasius friction factor is 
reasonable. 
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3 • A universal velocity profile is inadequate for use 
in the prediction of heat transfer characteristics 
of low Reynolds number turbulent flow. 

k. The use of the present low Reynolds number formula- 
tion leads to an accurate prediction of the Nusselt 
number for the thermal entrance find downstream, 
fully developed regions. For a constant wall heat 
flux, the results are described by the eigenvalue 
table and are correlated by the expressions 
presented in Chapter VI. 

5. The eigenvalue results (or the correlations for the 
eigenvalues and associated constants) provide the 
necessary information needed for the solution of 
the variable wall heat flux problem. 

6. It is not necessary to modify Reynolds analogy, 

t H - to obtain good heat transfer predictions 

for the conditions in this study. 
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APPENDIX A 

DETAILS OF EXPERIMENTAL VELOCITY STUDY 

Introduction 

This appendix presents details of the adiabatic 
velocity measurements taken at Texas Technological College, 
Lubbock, Texas, the data reduction procedure, and specifics 
of the uncertainty analysis . The uncertainty analysis 
includes consideration of the effect of neglecting probe 
corrections in relation to the uncertainty in the velocity 
and transverse position. The estimates are based on the 
expressions suggested by Daily and Hardison (27) in their 
impact probe literature survey. A brief discussion of 
these corrections is also included. 

Experimental Equipment 

Details of the experimental facility are presented 
in reference (25)* The major characteristics, plus modifi- 
cations performed for the present study, are documented in 
the main body of this work. The experimental facility is 
shown schematically in Figure A-l.' The probe installation, 
including the probe housing and wall static tap, is shown in 
Figure A-2. The radiation shield and insulation described 
in reference 25 were not in place when the present measure- 
ments were taken. 
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STATIC PRESSURE TAP 
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The manifold used in the pressure measurements was 
constructed from 1/k inch diameter Polypenco tubing and 
1/8 inch Swagelok tees. Six 1/8 inch Khitey and four 
1/8 inch Robbins valves with Swagelok connections were 
used. The six Whitey valves were connected to the tube 
wall static pressure taps. The remaining twenty-one wall 
taps were sealed by plugging the open end of the Swagelok 
connectors on the 1/8 inch wall pressure tap extensions. 
One-eighth inch diameter copper tubing, sealed at one end, 
was used for this purpose. The impact probe and its wall 
static tap were connected to the manifold through two 
Robbins valves. Approximately six to eight feet sections 
of 1/8 inch Polypenco tubing were used for these connec- 
tions . The high pressure and low pressure sides of the 
manifold were connected through two Robbins valves to a 
Flow Corporation Model MM3 Micromanometer with 1/k inch 
Tygon tubing. They were also connected directly' with 1/8 
inch Polypenco tubing to a Decker Model 306-2 differential 
pressure transducer and indicator. This arrangement 
allowed the pressure drop between any of three pairs of 
wall static taps, or the impact probe and its wall static 
tap, to be read by opening the appropriate valves on the 
Micromanometer and/or pressure transducer. The Micro- 
manometer has an accuracy of +0.0002 inches of butyl 
alcohol with a resolution of +0.0001 inch. The manifold 
was mounted on a 1/k inch sheet of plywood reinforced with 


■ y ■ -w . ** ' **-' '*'~t y*%^. 


h 'mww******' 


77 

slotted angle iron. It was clamped to a la rye sturdy 
table located next to the test section. The Micromanometer 
and transducer were also located on the table. 

A Flow Corporation Model HWB 3 Constant Current Hot 
Wire Anemometer was used to obtain qualitative fluctuating 
velocity component measurements . The output was observed 
on a Hewlett Packard 130-C oscilloscope , and a Sandborn 
4300 Series high speed recording oscilloscope was used to 
record u' traces. Due to an impedance matching problem 
between the high impedance hot wire equipment and the 
moderate impedance Sandborn, the Sandborn would not record 
the u’ output directly from the hot wire. Impedance 
matching was obtained by inserting a 20 microfarad capacitor 
in the positive leg of the hot wire output. Several resist- 
ances ranging from 100k$"^ to !.5mQwere also tried) but 
traces could not be recorded until the 20 microfarad 
capacitor was inserted in the system. The traces recorded 
in this manner were found to be superimposed on a sine wave 
of small amplitude and large period. Consequently , the 
traces obtained are only useful ns an indication of the 
oscilloscope output. 

The Flow Corporation hot wire probe, used to measure 
the u' fluctuations, was held with its axis coincident with 

the test sections by two circular oak spacers attached to 

the shaft of the hot wire probe. The spacers were 1.6 inch 

diameter wheels with three spokes and a 1/2 inch thick rim. 
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The spacers contained ample flow passages to prevent flow 
blockage. Access to the test section was gained through 
the exit of the 1.6l ID tube. The probe was positioned at 
the two stutions at which 11 * measurements were taken by 
means of a 1/k inch dowel attached to the Amphenol connector 
at the end of the probe. Scribed reference lines on the 
dowel and tube were used to reproduce the same relative 
alignment of the hot wire at the two stations. This 
apparatus was not present in the tube when veloc ity 
traverse datu were taken. 

Air was supplied by a compressor, capable of 
pumpxng ko cubic feet per minute, with the primary pressure 
varying from 85 to 125 psig. Flow control was maintained 
by a high capacity regulator and two Minneapolis Honeywell 
pressure regulators. The last regulator, a model 356529* 
was capable of maintaining control with 0.01 psig sensi- 
tivity and 0.013 psi drift in fifteen hours for pressures 
in the range of 0 to 25 psig. For low flow rates, all 

4 

regulators were used in series. Capacity limitations 
prevented the use of the Honeywell regulators for Reynolds 
numbers of 7000 and above, in which case they were by- 
passed . 

Bulk flow rates were measured with a Model 50 MC 
2-2P Meriam Laminar Flow Element. The reading was obtained 
directly in scfm on an inclined Meriam Manometer supplied 
with the unit. It was then corrected for pressure and 
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temperature effects* The combination has a stated accuracy 
of 1/2 of one per cent of reading over the entire 0.2 to 20 
scfm range* The temperature immediately downstream of the 
laminar flow element was measured with a Chromel -Alumel 
thermocouple. Pressure level was measured with either an 
8 inch inclined or GO inch vertical water manometer. 

A filtering system consisting of a Cuno cartridge 
filter, a silica gel dryer, and a Poro-Stone air filter was 
an integral part of the existing facility. To check the 
moisture content of the air, dew point measurements were 
obtained with an Electro Dryer dew point apparatus attached 
to the last pressure tap on the test section. The instru- 
ment yields a dew point temperature accurate to within £3 
degrees centigrade. 

The wall static pressure taps consisted of l/l6 
inch diameter holes in the tube, and 1/8 inch Inconel tube 
extensions 11-1/2 inches in length. The taps were con- 
structed by machining two concentric holes in the tube 
wall; a 1/2 inch milled hole approximately J/k of the way 
through the tube, and a l/l6 inch hole drilled through the 
remainder of the wall. The burrs were cut from the inside 
of the test section and the pipe was honed. Circular slugs, 
1/8 inch in thickness, to which the 1/8 inch Inconel 
extensions had been welded, were fitted into the 1/2 inch 
holes and welded to the test section. 




wwwb^^ .>*****#««*■ 


80 


impact Probe Calibration 

The impact probe was calibrate*! in an undisturbed 
laminar flow stream to determine the viscous effects dis- 
cussed in the literature survey of Daily and Hardison ( 27 ) • 
The probe error was determined by comparison of centerline 
measurements with the centerline velocity calculated from 
the measured bulk flow rate. The calibration, supplied by 
Texas Technological College, is presented in Figure A-3« A 
maximum overprediction of 0.02 per cent of the velocity is 
observed for local probe Reynolds numbers greater than 20 . 


Experimental Procedure 

Preliminary Check of Equipment and System 

Before data were collected, the experimental appara- 
tus was leak checked and preliminary pressure drop measure- 
ments were obtained to check the instrumentation. The test 
section exit was sealed and the system leak checked under 
50 psig pressure with M Snoop M leak detecting fluid. During 
this phase, all valves on the manifold were open, except 
those to the micromanometer , and the leads to the Decker 
transducer were disconnected and plugged. After a few 
initial leaks were discovered and sealed with U. S. Royal 
industrial adhesive (epoxy), no further leaks were found by 
this method. However, ths existence of minute leaks, 
estimated to be of the order of the accuracy of the 
micromanometer, was indicated by the initial pressure drop 



0.004 



Figure A-3* Probe Calibration 
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readings taken with flow. The readings obtained were 
highly erratic on the micromanomet er . In such cases, the 
micromanomcter would continuously drift in a given direc- 
tion. Tightening existing fittings and additional leak 
checking under pressure proved futile. Discussion with 
Dr. M. Davenport and Mr. L. Chance of Texas Technological 
College revealed, that while not mentioned in reference 
(50), II. Ferrell encountered the same problem. The 
problem was solved by coating all fittings and unions with 
paraffin. A systematic application of this technique to 
the present system eliminated the remaining leaks. 

Agreement between the Decker Pressure Transducer 
and the micromanomet er could not be obtained, for identical 
pressure drop measurements, in the subsequent instrument 
check out. The transducer proved to be erratic over short, 
unpredictable intervals, yielding both positive and nega- 
tive voltage readings. Instead of the stated fifteen 
minute warm up period, four to five hours sometimes elapsed 
before the zero balance could be obtained. In some cases, 
long periods of stability existed during which the output 
remained steady but differed from the micromanometer. 

These periods were terminated by instability of the 
instrument output and an inability to set the zero balance. 

An attempt to repair the transducer, following the 
instructions supplied with the instrument, proved futile. 

It was observed that, during the null balancing of the 
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micromanomot or , changes in the fluid level of the instru- 
ment induced oscillations between the niicromanometer and 
the pressure transducer. It was also observed that the 
opening and closing of valves on the manifold caused 
oscillations of the same nature. The transducer output 
was more sensitive to these oscillations than the micro- 
manometer . 

Based on the results of the preliminary adiabatic 
friction factor measurements, only the micromanometer was 
used to collect the pressure drop data for this investiga- 
tion. A separate set of data was also collected with the 
Decker transducer during one of its periods of relative 
stability. These data were not used in this work due to 
the following reasons: (l) disagreement with micromanometer 

results; ( 2 )- an inability to obtain a consistent calibra- 
tion of the transducer following the procedure originally 
used by Texas Technological College; and (3) the observed 
periods of instability discussed above. The data were 
recorded in data log book number 3 » on file with the 
Energy, Mass and Momentum Transfer Laboratory of the 
Department of Aerospace and Mechanical Engineering, 
University of Arizona. 

Fluctuating Axial Centerline Velocity Component > u» 

The qualitative turbulence measurements consisted 
of visual observations and chart recordings of the 
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centerline fluctuating axial velocity component, u', and the 
measurement of flow rates* All hot wire settings, the zero 
velocity current, penk-to-peak values of the calibration 
square wave and chart speeds were also recorded but were 
not used in the data reduction. Data were collected at two 
different axial positions for Reynolds numbers of 2000, 

2300, 2500, 2700, and 3000. These data were used to 
establish the limiting Reynolds number for which the flow 
at the probe was considered fully turbulent. Visual 
observation of the hot wire output is sufficient to indi- 
cate the Reynolds number at which the intermittency factor 
(17) has a value of one. Moreover, as mentioned previously, 
superposition of the u' traces on a sine wave of small 
amplitude and large period preclude quantitative evaluation 
based on the chart recordings. 

The hot wire data were obtained as follows. The 
hot wire probe and extension arm described in the equipment 
section were inserted into the test section exit and 
located at the probe station. The hot wire calibration 
square wave was then reproduced on the oscilloscope and 
5 ■mdborn recorder to assure that repr esent at -on of the hot 
wire output was being obtained. Following this preliminary 
step, a desired flow rate was set and the gas temperature 
and pressure recorded, as described in the next section. 

The hoi wire output was then observed for several minutes 
on the oscilloscope and the observed flow patterns and 
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axial positions were indicated in the data book* The hot 
wire output was then fed into the Sandborn Recorder and a 
trace of the output obtained on light sensitive recording 
paper* Without changing the flow settings, the probe was 
then moved to the second axial position and the observation 
and recording process repeated. Next, a new flow setting 
was established and the entire procedure repeated for the 
two axial positions* 

At the completion of the test, the flow was turned 
off by closing the high capacity flow regulator. The hot 
wire calibration square wave was then reproduced again on 
the oscilloscope and Sandborn recorder to assure that the 
output of the hot wire equipment was still being obtained. 

Adiabatic Flow Runs 

The adiabatic flow runs consisted of measurements 
of wall friction factors for three different pressure tap 
pairs, impact (velocity) probe traverse data, and flow 
rates. Data were collected at bulk Reynolds numbers of 
3020, 4080, 5010, 7030, 10,100, and 15,000. All runs were 
performed during the evening hours from 7 p.m. to 3 a.m. 
This procedure was necessitated by the extreme sensitivity 
of the micromanometer, which was affected by the rapid 
changes in atmosphere conditions that occurred during the 
day, as well as by air currents induced by the evaporative 
cooling system. 
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The procedure employed for a given flow run was as 
follows. Approximately two hours before testing began, the 
evaporative cooling system was shut off. At the end of the 
two hour period, the air compressor was turned on with all 
three of the flow (pressure) regulators, and the main flow 
valve closed* When the pressure tanks were filled with 
compressed air at 125 psig, the regulators were adjusted to 
give the approximate desired flow rate. The barometric 
pressure was measured. Based on temperature and pressure 
measurements of the flowing gas, the laminar flow element 
output necessary to produce the desired Reynolds number 
was computed. This output was then set on the laminar flow 
element by fine adjustment of the regulators. For low flow 
rates (Reynolds numbers below 7000) all three regulators 
were used in series. For high flow rates, the two low 
capacity regulators were bypassed. 

The time was then recorded and the dew point 
temperature of the air determined. The level of the 
micromanometer was checked with the integral circular 
spirit level and adjusted if necessary. The micro- 
manometer zero (AP = 0) was then determined. Adiabatic 
wall friction data were then collected for the three 
different wall pressure tap pairs. These data consisted 
of the pressure drops obtained by null balance of the 
micromanometer, in inches of butyl alcohol, and the alcohol 
temperature. For the first and third pressure tap pairs 
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the laminar flow element output in scfm, the gas tempera- 
ture and pressure, and the barometric pressure and 
temperature were also recorded. 

The velocity traverse was then initiated. Follow- 
ing a check of the micromanometer zero, the difference 
between the probe impact pressure and wall static pressure, 
in inches of butyl alcohol , and the temperature of the 
alcohol were recorded for a given micrometer reading of 
the traversing mechanism. A tv averse, consisting of 17 to 
20 positions, was performed with the probe moving from the 
wall towards the centerline. Small, radial changes were 
employed near the wall, increasing in size from approxi- 
mately 0.01 inches to 0.075 inches in the vicinity of the 
centerline. Starting with the first probe setting, the 
laminar flow element output in scfm, gas temperature and 
pressure, and the barometric pressure and temperature were 
recorded for approximately every second or fifth probe 
setting. Due to the high sensitivity of the micromanometer 
and its null technique of measurement, the flow runs were 
quite lengthy, averaging from 2-1/2 to 3-1/2 hours per run. 

At the conclusion of a traverse, the micromanometer 
zero was rechecked and the time recorded. If time per- 
mitted, a new flow setting was set and the process 
repeated. Otherwise, the system was shut down. The shut 
down sequence consisted of placing the micromanometer and 
laminar flow element off line and opening their bypass 
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valves. The main flow valve and the regulators in use 
were then closed, the compressor turned off, and the 
evaporative cooling system turned back on. 

Data Reduction 

Fluctuating Axial Centerline Velocity Component, u* 

For the hot wire anemometer measurements of the 
fluctuating axial centerline velocity component, data 
reduction was not attempted. The traces were examined 
visually on the oscilloscope to yield a qualitative measure- 
ment of the Reynolds number above which fully turbulent flow 
existed at the measurement station. Representative traces 
at the probe location, recorded by the Sandborn recorder, 
are presented in Figure A-^l. The decay of the laminar 
segmert observed on the middle trace is probably due to 
the electronic circuitry. On the Hewlett Packard Oscil- 
loscope, laminar flow appeared as a straight, uniform, 
smooth line. 

Adiabatic Flow Runs 

Values for the viscos ity, |i, and the compres- 
sibility factor, Z, used in the velocity and friction data 
reduction were taken from Iiilsenrath et al . (51). Two way 
interpolation on pressure and temperature was employed in 
all cases. Measurement of extremely low dew point 


R* ■ 2,000 
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Figure A-4. Representative u' Traces 
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temperatures, in the vicinity of -*» B degrees Fahrenheit tv r 
all tests, justified the use of dry air properties. 

The mass flow rate was calculated from the scfm 
output of the Meriam Laminar Flow Element , with corrections 
for temperature and pressure effects. The density was 
calculated from the equation of state , including the 
compressibility factor. For the Reynolds number, the 
viscosity was evaluated at the gas temperature and pressure. 
The wall shear stress was obtained by a force balance 
between any two pressure taps with the pressure drop 
measured by the micromanometer. The friction factor was 
evaluated from the Fanning friction factor expression. 

The velocity was calculated from the expression 


u(y ) 



2g c AP 

~P~ 


( A-l ) 


The pressure drop between the impact probe and its wall 
static tap, A P, was obtained from micromanometer readings 
in inches of butyl alcohol* The alcohol density was 
evaluated as a function of temperature from the information 
supplied with the micromanometer. An estimation of the 
effect of neglecting the contribution of the fluctuating 
velocity components in the above equation is included in 
the uncertainty analysis, presented in the next section. 

The nondimensional velocity, u , and distance, y , are 


u* . ,flM. 

V r w*c *p 
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♦ yVr « /p 
, y = — ^ 

An forming thoso two groups, the properties w ere evaluated 
at the measured gas temperature and pressure. The experi- 
mental values of wall shear stress were used. 

The local probe Reynolds number was calculated for 
each probe position from the expression 


Re 


P 


u ( r ) D 
E 

v 


(A-3) 


Corrections from the probe calibration presented in Figure 
A-2 were unnecessary since the probe Reynolds number always 
exceeded twenty. 

The bulk velocity was determined in two different 
ways: directly from measurements of the bulk flow rate, and 

from a graphical integration of the velocity profile, u(y) . 
A comparison of the t?m values is presented in Appendix C. 


Uncertainty Analysis 

Important considerations of the uncertainty 
analysis, which is an estimate of the cumulative effect of 
individual errors, are presented below. The method employed 
was that presented by Kline and McClintock (30). The 
uncertainties in the directly measured quantities are 
presented in Table A-l. These values were based on con- 
sideration of the manufacturer's specifications (where 
available), the estimated reading error, and experience. 
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Table A-l. Velocity Study: Uncertainty Intervals of 
Measured Quantities 


Item 


Uncertainty 


Luminur Flow Element (Flow Rate) 

Micromanometer (Pressure Drop) 

Gas Pressure (above barometric) 

Barometric Pressure 

Probe Location (Micrometer) 

Temperature 

Tube Diameter 

Distance Between Pressure Taps 
Viscosity 


+ 1 / 2 % of reading 
♦0.00002 psi 
♦ 0.005 in » 2 0 
♦0 .0204 psi 
♦0.0005 inch 
♦^•F 

♦p.001 inch 
♦ 0.006 inch 
♦0.01 lbrn/ift Bile 
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The uncertainties in the deduced results are presented with 
the tabulated results in Appendix B. 


Probe "Error” Uncertainties 

An estimate of the magnitude of the error caused by 
neglecting probe corrections was made. As discussed in 
Daily and Hardison (27) * the effect of turbulence on impact 
probe readings* referred to a datum pressure measured at the 
wall* is expressed by the relation 





(A-4) 


The velocity is seen to be 


where 



a «c (p - V 
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+ 
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A common procedure, employed by Laufer (52) and 
others who have made probe corrections* has been to 
neglect the turbulence effect on the wall static probe (the 
last two terms) in using the above correction. Complete 
neglect of correction C is equivalent to introducing a 
maximum error * W * 


in u 
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or, conversely, a maximum percentage error in u of 
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= V l + 2 *c AP " 
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This error is now treated as an effective uncertainty in 
the velocity. The total uncertainty in the velocity 
becomes 


W , W A , 2 . W„ 

u / 1_ AP v /IP 

u " L 2 AP ' v 2 p 
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(A-9) 


W 

To estimate the importance of -p, a logarithmic 

extrapolation of C was performed. Sandborn's ( 26 ) 

normalized turbulence measurements, obtained at Reynolds 

numbers of 25,000 and 50,000, were used in this evaluation 

W 

Q 

of che correction term. In Table A-2 values of — — are 
compared with the sum of the other three uncertainty terms 
for a Reynolds number of 4080. The other terms dominate. 
Subsequently, this term was neglected in estimating the 
uncertainty of the velocity. 


Effective Center Displacement Effect 

Insertion of a finite probe into a shear flow, or 
flow near a boundary, causes an "'effective center dis- 
placement" of the probe. In effect, the probe may be 
considered to measure the stagnation pressure of a 
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Table A-2. Importance of Turbulence Velocity Correction 


y/r 
J w 

♦ 

y 

Re » 4080 
V 2 

Term (*~) 

All other terms , 
equation A-9 

0.0472 

.?44 

1.42 

11.91 

0.1093 

15.61 

0.307 

6.64 

0.2335 

33.36 

0.0984 

3.61 
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streamline which is shifted from its upstream course due to 
the presence of the probe. Normally the displacement of 
the streamline has been considered as away from the wall. 
Thus, the velocity measured by the probe is closer to the 
wall than the location of the probe centerline. However, 
streamline displacements toward the wall have been 
reported (27) . No comprehensive data are available upon 
which to base a correction (27) • To estimate the possible 
error caused by this effect, Fage (29) analytically extended 
experimental results obtained by Young and Maas (see 27) • 

In the present work the effective center displacement , A » 
given by Fage’s expressions, 

d 2 

A* 90(0.131 ♦ 0.o8 fc * d./d )y (-“) 

x o r 

w 0 «e y/r w « 0.15 

and (A-10) 

d 2 

A « 16(0.131 + 0.082 d./d Hr - y)(-~) 

1 ° W r w 0.15 S y/r wS l 

was considered as one of the errors in the probe position. 
Thus, the total uncertainty in the position of the velocity 
was evaluated from 

W _ /~\j 2 T 

y * y y 

where W is the Uncertainty interval of the probe micro- 
meter positioning mechanism. 


97 


Effect of Slot in Tube Wall 

In the vicinity of the wall a large discrepancy is 
observed between the velocity data and accepted viscous 
profiles. The uncertainty analysis is completely unable to 
account for this divergence, which is believed to be caused 
by the small rectangular slot milled in the tube wall to 
accommodate the probe traversing mechanism. The width of 
the slot is approximately 1 per cent of the tube circum- 
ference. It could cause changes in the turbulence level 
and/or an effect comparable to the effective center 
displacement . 

Haugen and Dhanak ( 31 ) have presented experimental 
velocity profiles for boundary layer flow over rectangular 
slots, at estimated boundary layer Reynolds numbers of 
12,900, 38,600, and 64,000. An attempt was made to predict 
the y + values at which disagreement between these profiles 
(presented ii. graphical form) and Spalding's law of the 
wall profile (53) would first occur. Agreement was found 
for the smallest distance from the wall which could be 
accurately plotted on the graphs. Unfortunately, the 
distance represented a sizeable y value (approximately 
60), and a comparison near the wall was not possible. 
Consequently, the results only indicate a point at which 
the effect of the slot is known to be negligible. The 
ratio of this distance over the boundary layer thickness 
was plotted against the boundary layer Reynolds number, for 
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each of the three values above, on logarithmic coordinates. 
Extrapolation to the Reynolds numbers of the present study 
yielded a value of y + of approximately forty for all of the 
tests . 

From a qualitative examination of the present data, 
it appears the slot effect is no longer felt beyond y/r 

W 

ranging from approximately 0.09 to 0.05 (y + * 9»7 to 21.8), 
corresponding to the lowest and highest flow rates of the 
present study. This range of y/ r w represents the locations 
at which agreement between the velocity data and accepted 
viscous profiles (where appropriate) and/or strong changes 
in slope, du /dy , are first observed. As indicated 
previously, the effect can be conservatively considered 
negligible by y + * ko . 

Errors Due to Selection of Wall Shear Stress 

The values calculated for y + and u + are directly 
dependent on the expression or method utilized to evaluate 
the wall shear stress since it appears in the shear velocity 

( A-12) 

used in their definitions. Two approaches have been taken 

to this problem: (a) use of an analytical expression for 

\he friction factor, or (b) evaluation of T from experi- 

w 

mental data taken during the velocity investigation. 

Senecal (4) employed a friction factor evaluated from 



Blaaius* expression. The second alternative was used in 
the present study, but, when reduced by both methods, the 
resulting values of u* show a difference of only 0.5 per 
cent. At low Reynolds numbers this difference may exist 
in comparisons between the present formulation and the 
McEligot, Ormand and Perkins profile (9)t which was based 


on Senecal's data 
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APPENDIX D 

EXPERIMENTAL EDDY DIFFUSIVITY VALUES 


Values of the eddy diffusivity determined from the 
present velocity profile data (via equation D-l below) are 
presented in Figure B-l . These profiles were employed as 
an aid in the selection of the values of y* and % which 
appear in equation 3 “ 5 » The method of formulation and 
subsequent use of these values are discussed in the 
section, Semi-empirical Profile. 

Eddy diffusivity distributions predicted from the 
combined form of Reichardt's wall and center laws (2*l), 
with his original values for y* and are presented for 
Reynolds numbers of 3020 and 15*000 in the figure. 

It is observed that the experimental values tend to 
fall towards negative values for values of y/r in excess 
of 0.6. This drop is a direct consequence of the inability 
to obtain accurate experimental results for the eddy 
diffusivity from the defining equation 



(B-l) 


m 


which was used to obtain from the velocity data. 


Analytically* this expression tends to an indeterminate 
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Figure B-l • Experimental Eddy Diffusivity Results 
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form as the centerline is approached. However* in data 

f + 

reduction the percentage uncertaint> in du /dy becomes 
largo despite precise velocity determination a aero value 
of the slope is seldom obtained. Accurate experimental 
values of the eddy diffusivity in the vicinity of the 
centerline have been determined from hot wire turbulence 


measurements 
finite value 


by several investigators 


of ■ at the centerline 

-u>i 


(32). These yield a 
of approximately 0.07 


rather than the steep drop observed in Figure 13-1. 


The values shown for less than 0.1 are probably 

w 

of little use due to the errors in the velocity near the 


wall t as discussed earlier. Conservatively! the most valid 

data thus fall in the region 0.2 « S: 0.7* 

w 


APPENDIX C 


ADIABATIC BULK VELOCITY COMPARISON 


As an overall check on tho experimental velocity 
profiles , the bulk velocity was determined by graphical 
integration of the velocity data. The integral appearing 
in the definition of the bulk velocity was determined by 
plotting the profile and evaluating the area with a 
plnnimeter. In Table C-l the results are compared with 
values of bulk velocity determined from the measured mass 
flow rate. The approximate "percentage over estimation" of 
the bulk velocity, obtained in the integration, due to the 
effective center displacement effect is also included. 

This estimation was obtained using Page's expression (29) , 


- V. 


correct 

v b 


g r d 2 

- 190 + 190 (-£) ♦ 

P v b 


120 


w 


d . d 
X o 


w 


(C-l) 


However, Fage's formulation was not based on comprehensive 
data (27) and must be considered only as an indication of 
the error in bulk velocity obtained by integration. More- 
over, since both positive and negative effective center 
displacements have been reported, this approximation 
probably should be considered as representing a possible 
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Table C-l. Comparison of Bulk Velocity Calculations 


Re 

V. . 
bx 

integration 

ft/sec 

V bm 

measured 

ft/sec 

% difference 
v bi. - V b„ 
v b± 

% M overestimation M 
in V| } integration 
(Fage) 

3,020 

4.53 

*1.37 

+ 3.53 

+ 1.**1 

*k,o 8 o 

5.35 

5.33 

+ 0 . 3*12 

±1 .3*1 

5,010 

6.97 

7.03 

- 0.86 

+.1 .305 

7,030 

9.55 

9.65 

-1 . 0*1 

+ 1 . 2 *i 

10,100 

13.99 

1*1.19 

-1 .*13 

♦ 1 .16 

13,000 

19.99 

20 . kl 

- 2.10 

* + 
• 

O 

0= 
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bound on the error in the bulk velocity without regard to 
direction . 


APPENDIX D 


TABULATED VELOCITY DATA AND UNCERTAINTIES 
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Adiabatic Velocity Data: Air 

Hun 1 He s 4077 

f * 0.00982, 0.00929, 0.01052 

x/D range 234-251, 247-264, 259-275 

T a 88 # F , P a 27.015 in Ilg. abs. 
av 7 av ■ * 


y 

u 

X 

X 

<6y+)/y* 

<6'> + )A> + 

inches 

ft/sec 

y 

u + % 

> uncertainty °/ 

'0 uncertainty 

0.004 

2.026 

0.710 

4.957 

13.0 

34.6 

0.012 

2.219 

2.130 

5.430 

5 * 5 

28.9 

0.028 

2.718 

4.969 

6.651 

4.1 

190 

0.038 

3.450 

.6.744 

8.442 

3.9 

12.2 

0.044 

3.624 

7.809 

8.868 

3.8 

11.1 

0.053 

3.680 

9.406 

9.005 

3.8 

10.7 

O.063 

3.949 

11.181 

9.663 

3.7 

9.4 

0.073 

4.201 

12.956 

10.279 

3.7 

8.4 

0.088 

4.620 

15.618 

11.304 

3.7 

7.1 

0.129 

5.695 

22.894 

13.933 

3.6 

5.0 

0.148 

5.907 

26 . 266 

14.453 

3.3 

4.7 

0.188 

6.277 

33.365 

15.359 

3.0 

4.3 

0*238 

6.502 

42.239 

15.909 

2.9 

4.1 

0.313 

6.841 

55.550 

16.737 

2.7 

3.9 

0.391 

7.048 

69.393 

17.244 

2.7 

3.7 

0.471 

7.361 

83.591 

18.010 

2 . 7 

3.6 

0.563 

7.581 

99.918 

18.548 

2.6 

3.5 

O.638 

7.635 

113.229 

18.6 80 

2.6 

3.4 

0.713 

7.662 

126.539 

18.746, 

2.6 

3.4 

0.788 

7.688 

139.850 

18.811 

2 . 6 

3.4 
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y 

inches 

0.013 

0.028 

0.038 

0.044 

0.053 

O.O63 

0.073 

0.088 

0.113 

0.138 

0.163 

0.188 

0.238 

0.288 

0.338 

0.388 

0.463 

0.563 

O.638 

0.713 




Adiabatic Velocity Data: Air 

Run 2 Re * 3017 
f = 0.01216, 0.01056, 0.01031 


x/D 

range = 

234-25? , 

247-264, 259 

-275 

T 

av 

= 91. 4°f, 

P * 26.970 in Hg. 

av ' w 

nbs. 

u 

a, 


(6> + )/y + 

« 5 « + ) / u * 

ft/sec 

▼ 

y 

T 

u 

% uncertainty 

% uncertainty 

2.131 

1.770 

6.712 

6.2 

31.7 

2.649 

3.813 

8.344 

3.2 

20.7 

2.726 

5.175 

8.386 

5.1 

19.6 

2.801 

. 5.992 

8.821 

5.0 

18.7 

2.944 

7.217 

9.274 

5.0 

17.0 

3.014 

8.579 

9.49?. 

4.9 

16.2 

3.276 

9.941 

10.319 

4.9 

13.9 

3.577 

11.984 

11.268 

4.9 

11.9 

3.908 

15.388 

12.310 

4.9 

10.2 

4.310 

18.793 

13.576 

4.7 

8.7 

4.678 

22.197 

1^.733 

4.5 

7.7 

4.851 

25.602 

15.279 

4.4 

7.3 

5.l4o 

32.410 

16.190 

4.3 

6.7 

5.376 

39.219 

16.932 

4.2 

6.4 

5.636 

*;6 .028 

17.758 

4.2 

6 .0 

5.711 

5i.837 

17.988 

4.2 

5.9 

5.854 

63.051 

18.437 

4.2 

5.8 

5.924 

76.668 

18.658 

4.2 

5.7 

6.028 

86.882 

18.985 

4.2 

5.6 

6.130 

97.095 

19.305 

4.2 

5.5 

6.163 

107.309 

19.411 

4.2 

5.5 


0.788 






y 

inches 

0.013 

0.020 

0.028 

0.033 

0.038 

0.044 

0.033 

0.063 

0.073 

0.088 

0.113 

0.138 

0.188 
0.238 
0.288 
o . 388 

0.463 

0.563 

0.638 

0.713 


W0P 
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Adiabatic Velocity Data: Air 

Run 3 he * 10070 


f = 0.00780, O.OO783, 0.00774 
x/D range * 234-251, 247-264, 259-273 


T 

av 

= 87. 7°F , 

P av " 

. 400 in Ilg. abs 

• 

u 

4 . 

4 . 

<6y*)/y* 

<6"*>/** 

ft/sec 

V 

11 °/o 

uncertainty % 

uncertainty 

1.680 

5.087 

1.893 

4.8 

49.5 

5.538 

7.827 

6.238 

3.8 

4.6 

6.812 

10.957 

7.673 

3-3. 

3.1 

8.085 

12.914 

9.107 

3.2 

2.2 

8.307 

14.871 

9056 

3.1 

2.1 

9.653 

17.219 

10.875 

3.0 

1.6 

10.144 

20.741 

11.426 

3.0 

1.5 

IO.892 

24.654 

12.268 

2.9 

1.3 

12.002 

28.567 

13.519 

2.9 

1.1 

12.399 

34.437 

13.966 

2.9 

1.0 

12.972 

44.221 

14.611 

2.8 

1.0 

13.365 

54.004 

15.279 

2.5 

0.9 

14.261 

73.571 

16.063 

1.9 

0.9 

14.679 

93.138 

16.534 

1 .6 

0.8 

15.232 

112.704 

17.157 

1.4 

0.8 

16.307 

151.838 

18.368 

1.3 

0.7 

16.722 

181.188 

18.836 

1.2 

0.7 

17.292 

220.321 

19.477 

1.2 

0.7 

17.821 

249.671 

20.073 

1.2 

0.7 

17.933 

279.021 

20 . 200 

1.2 

0.7 

17.990 

308.371 

20.263 

1.2 

0.7 


0.788 
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Adiabatic Velocity Data: Air 

Hun 4 Re = 5010 
f = 0.00927, 0.00964, erratic 
x/D range = 234-251, 247-264, 259-275 
T av = 8l * 6 3 # P, P ftV = 26.979 in Hg. abs. 


y 

u 

i, 


(<5> + )/y + 

«5u + )/« + 

inches 

ft/sec 

y 

U /'i 

s uncertainty % 

uncertainty 

0.013 

1.105 

2.754 

2.310 

5.0 

114.9 

0.038 

2.707 

8.051 

5.658 

3.5 

19.2 

0.044 

3.552 

9.322 

7.425 

3.5 

11.3 

0.053 

3.933 

11.229 

8.221 

3.4 

9.2 

0.063 

4.420 

13.348 

9.239 

3-4 

7.4 

0.073 

5.525 

15.467 

11.549 

3.3 

4.9 

0.088 

5.882 

18.645 

12.295 

3.3 

4.4 

0.113 

6.444 

23.942 

13.468 

3.3 

3*8 

O.138 

6.931 

29*239 

14.486 

3.0 

3.4 

0.188 

7.386 

39.832 

15.437 

2.5 

3.1 

0.238 

7.709 

50.426 

16.113 

2.3 

2.9 

0.288 

7.943 

61.020 

16.603 

2.2 

2.8 

0.388 

8.319 

82.207 

17.387 

2.1 

2.7 

0.463 

8.560 

98.098 

17.891 

2.1 

2.6 

0.563 

8.841 

119.286 

18.478 

2.0 

2.5 

O.638 

8.955 

135.176 

18.717 

2.0 

2.5 

0.713 

9.051 

151.067 

19.000 

2.0 

2.4 

0.788 

9.135 

166.957 

19-094 

2.0 

2.4 
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Adiabatic Velocity Data: Air 

Hun 3 Re * 7025 
f « 0.00849, 0.00853, 0.00843 



x/D 

range = 2 

3<t-25i, 

247-264, 259- 

275 


T 

av 

= 79.75 

* **av 

27.107 in Hg. 

abs. 

y 

u 

+ 

y 


((5y + )/y* 

«5u + )/u + 

inches 

ft/sec 

¥ 

u 

% uncertainty 

% uncertainty 

0.013 

1.270 

3.708 

1.994 

4.8 

86.2 

0.038 

4.538 

10.840 

7.119 

3.2 

6.8 

0.044 

3.927 

12.551 

9.298 

3.1 

4.1 

0.053 

6.226 

15.119 

9.767 

3.1 

3.7 

0.063 

6.989 

17.971 

IO.966 

3.0 

3.0 

0.073 

7.860 

20.824 

12.331 

3.0 

2.5 

0.88 

8.236 

25.103 

12.921 

3.0 

2.3 

0.113 

8.896 

32.234 

13.956 

2.9 

2.0 

0.1fl8 

9.339 

39.366 

14.651 

2.6 

1.9 

O.188 

9.986 

53.629 

15.667 

2.1 

1.7 

0.238 

10.285 

67.892 

16.136 

1.8 

1.6 

0.288 

lo.84o 

82.155 

17.006 

1.7 

1.5 

0.388 

11.242 

110.681 

17.637 

1.5 

1.5 

0.463 

11.613 

132.075 

18.219 

1.5 

l .4 

0.563 

12.106 

160.601 

18.993 

1.4 

1.4 

0.638 

12.239 

181.996 

19.201 

1.4 

1.4 

0.713 

12.419 

203.390 

19*484 

1.4 

1.3 

0.788 

12.565 

224.784 

19.712 

1.4 

1.3 
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Adiabatic Velocity Data: Air 

Run 6 Re « 15031 
f s 0.00706, 0.0070(1, erratic 
x/D range * 234-251, 247-264, 259-275 


T av s 78.50T, P av = 27.57(1 in llg, abs. 


y 

u 

▲ 

▲ 

<6y*)/y* 

<6»*)/u* 

inches 

ft/sec 

▼ 

y 

u* % 

uncertainty % 

uncertainty 

0.013 

4.361 

7.210 

3.595 

4.7 

7.2 

0.018 

6.48o 

9.983 

5* 3^2 

3.9 

3-3 

0.023 

8.902 

12.756 

7.338 

3.5 

1.7 

0.028 

10.645 

15.329 

8.775 

3.3 

1.2 

0.033 

11.876 

18.302 

9.790 

3.2 

1.0 

0.038 

13.143 

21.075 

10.834 

3.1 

0.8 

0.044 

14.962 

24.403 

12.334 

3.0 

0.7 

0.053 

15-483 

29.394 

12.763 

2.9 

0.6 

0.063 

16.317 

34.940 

13.451 

2.9 

0.6 

0.088 

17.948 

48.805 

14.795 

2.8 

0.5 

0.138 

19.237 

76.335 

15.858 

2.5 

0.5 

0.188 

20.309 

104.266 

16.741 

1.9 

0.4 

0.238 

21.131 

131.996 

17.419 

1.5 

0.4 

0.388 

23.102 

215.187 

19.044 

1.2 

0.4 

0.463 

23.795 

256.782 

19.615 

1.1 

0.4 

0.563 

24.878 

312.243 

20.507 

1.1 

0.4 

O.638 

25.241 

353.838 

20.807 

1.1 

0.4 

0.788 

25.770 

437.029 

21.242 

1.1 

0.4 


APPENDIX K 


DETAILS OF EXPERIMENTAL HEAT TRANSFER MEASUREMENTS 

Experimental Apparatus 

The experimental gas loop is shown schematically 
in Figure E-l . The gas supply consisted o*’ four largo 
commercial compressed gas cylinders. The test gas was 
reduced to operating pressures by the two pressure 
regulators and flowed through the primary and secondary 
flow meters to the inlet mixer where its bulk temperature 
was measured. It then passed through the vertical , 
resistively heated, circular test section and entered a 
shell and tube heat exchanger where the hot gas was cooled. 
It was exhausted through the downstream control valves to 
the atmosphere. Measurements of the test section pressure, 
overall pressure drop, wall temperature profile, and 
approximate exit bulk temperature were obtained. The test 
section was enclosed in a vacuum chamber to reduce the heat 
loss and allow localized heat loss calibration. The vacuum 
reduced the time necessary to reach equilibrium with 
heating . 
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ilk 



. Figure E-i. Heat Transfer Apparatus 
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Teat Section and Vacuum Chamber 
General Details 

The teat section was fabricated from a 0.250 OD 
Inconel tube with a 0.010 inch wall. The exit end was 
brazed to a 3/8 inch thick wall stainless steel tube 
approximately 10 inches in length. The entire section was 
hung from a vertical stainless steel, flange, by means of a 
Conax fitting, inside a 6 foot long, 6 inch diameter glass 
vacuum chamber. The glass cylinder was located on top of a 
large glass cross of equal diameter, the three free arms of 
which were sealed by stainless steel flanges. A flexible 
5/8 inch 0D high pressure hose was used to connect the 
.test section inlet to a fitting, welded to one of the 
flanges. The flexible hose allowed the test section to 
expand under its own weight when heated. Vertical align- 
ment of the test sectior was ensured by Teflon and 
plexiglass guides in unheated regions. 

Particulars of the test section are presented in 
Table E-l. Thermocouples of $6 gage (0.005 inch) Chromel- 
Alumel wire were spot welded to th A test section using a 
Weldmatic Model 1015-C Capacitor-Discharge spot welder. 
Twenty gage (0.0321 inch) thermocouple wire was used 
between the fine wires and the two isothermal junction 
boxes, where a conversion to copper wire was made. The 
isothermal box for thermocouples 1 through 13 and the inlet 
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Table E-l. Test Section Information* 


Material 



Inconel 600 

Inside diameter 


0.2302 

Outside diameter 


0.2502 

Hydrodynamic 

entrance length 


23.0 

Heated length 


25.0068 

Pressure tap 

holes 


o.oko 

Item No. 

Location 

Item No. 1 

Location 

Location of 

pressure taps 



1 

- 1 .6862 

2 

27.2068 

Wall Thermocouple locations 



1 

-3.3051 

9 

5.9620 

2 

-2.1798 

10 

7.9610 

3 

-0.5962 

11 

11.9612 

k 

0.0 (L. E.) 

12 

19.9498 

5 

0.2781 

13 

22.0068 

6 

0.9387 

Ik 

25.0068 (U. E.) 

7 

1.8718 

15 

25.3010 

8 

3.6996 

16 

26.6905 

(L. E. = 

Lower Electrode 

, U . E . as Upper 

Electrode ) 

Glass thermocouple locations 

with reference 

to lower 

electrode 




1 

30.2568 

k 

1.2568 

2 

22.0068 

5 

- 9.7432 

3 

12.2568 




*A11 dimensions in inches 
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bulk temperature was located inside the vacuum chamber. 

Its isothermal reference temperature was given by a 
precision mercury-in-glass thermometer. For all other 
thermocouples, the reference junction was an ice bath. 

Test Section Diameter 

The transverse dimensions of the test section were 
measured by several different means, using tubing obtained 
from the same mill run. Results are presented in Table E-2. 
Values chosen for the ID and CD were 0.2302 and 0.2502 
inches . 

Instrumentation and Equipment 

The supply gas pressure was reduced to the desired 
operating pressures by means of a two stage Victor regulator, 
Model VTS 431-250, and a Wilkerson Model 2001-4 pressure 
regulator, connected in series. The flow rate was measured 
by a Brooks Model 1112 A Variable area flowmeter (helium), 
or a Brooks Model 150 Sho-Rate flowmeter (air), with stated 
accuracy of +1 per cent of full scale reading. A Foxboro 
Model 13-A differential pressure transmitter, with integral 
orfices, served as a secondary check. 

The thermocouple signals were fed to a Honeywell 
Electronik 18 self balancing indicating potentiometer 
through an adjustable range unit. The Electronik 18 , 
calibrated by Honeywell, has an accuracy of +0.25 per cent 
of span (5mv) and the range unit has an accuracy of +0.1 
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Table E-2. Measurements of Tube Dimensions 




ID 

WALL 

OD 

Test 

Instrument s 

inches 

inches 

inches 


1 Starret micrometer & 
small hole gage 

(Average--3 readings) 0.2302 -- 0.2500? 


2 Millers Falls micrometer 

& small hole gage 0.230 

3 Craftsman #40l6l-DJ 

helio vernier calipers 0.228 


0.011 


0.250 


0.250 


k Pratt and Whitney 
electrolimit gage 


0 . 252 (max) 
0.2507 (min) 


5 Unitron Model MEC 
inverted stage 

metallurgical microscope 0.23043 0.010250 0.24888 


Tests 2, 3» and 4 used the same sample 
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per cent of scale suppression. The accuracy of thi® 
combination is greater than the .accuracy of the thermo- 
couples used. Premium grade Chromel -Alumel thermocouple 
wire, accurate to + 2°F below 530°^ and +P»375 per cent of 
reading above 530°*' * was used throughout the system. The 
test section thermocouples were spot welded to the outside 
of the test section; those located on the heated length 
were calibrated in place for thermocouple conduction error 
as described in the sec tion--Experimental Procedure. The 
entering gas bulk temperature was measured in an inlet gas 
mixer. A mercury-in-glass NBS Certified Bomb Calorimeter 
thermometer, with a least scale division of 0.05°F, was 
used to measure the temperature of* the isothermal junction 
box located inside the vacuum chamber. 

The pressure drop was measured either on a Meriam 
Model 3^FB2 Micromanometer, on a 10 inch inclined manometer, 
or on a 60 inch water manometer. These instruments were 
accurate to about 0.001, 0.01, and 0.05 inches, respec- 
tively. For the heating runs, which were at "high" 
pressures, the micromanometer was not employed and the test 
section and flow meter pressure levels were measured on a 
150 psi Heise gage, accurate to 0.15 psi. For some 
preliminary adiabatic runs either a water manometer or a 
mercury manometer was used to measure the pressure levels. 

A Welch Model 1215 mercurial barometer was used to measure 
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the atmospheric pressure, which was corrected for tempera- 
ture and local gravity effects. 

The test section voltage drop was measured with a 
Weston Model 3^1 Voltmeter, accurate to 0.25 per cent of 
full scale. (For the air data « Fluke Model 883A AC/DC 
Null Voltmeter, with uncertainty of +0 . 1 per cent of read- 
ing, was used.) The voltage measurements provided a check 
of the overall test section resistance and the local energy 
generation. The current was measured with a Weston Model 
370 Ammeter and Weston Model 327-Type 2 Current Trans- 
former. The ammeter and current transformer have a 
combined uncertainty of about 0.25 per cent of full scale. 
These voltmeters and the ammeter were used in the calibra- 
tion for the test section resistance. The calibration is 
discussed in Chapter VII. 

Power for the electrical heating was provided by a 
surplus 3&VA GE Model 6lC-75 transformer. The 110 V, 60 
cycle AC input of the transformer was stabilized to within 
+0.01 per cent by a Sorenson Model 250-1 voltage regulator. 
The input current was controlled by an Ohmite Model VT 20 NB 
Variable Transformer. Above approximately 80 amperes a 
Lincoln Model TM-500 AC/DC welder was used in place of the 
above power supply package. 

The vacuum pumps, which maintained the system below 

.5 

1.5 > 10 torr, consisted of a 2 inch Consolidated Vacuum 
Corporation Model PMC-115 diffusion pump and a Welch 
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Duo-Seal Model l / i02U mechanical pump. Vacuum pressures 
were measured with either a Pirnni gage tube and a 
Cojisolidat ed Vacuum Corporation Type GP-l4tO indicator, or a 
General Vacuum Industries Model 700 cold cathode vacuum 
gage . 


Experimental Procedure 

Preliminary Tests and Calibrations 

Before collection of the heat transfer data was 
undertaken, preliminary tests and calibrations) designed to 
increase the accuracy of the data and to assure correct 
operation of the equipment) were performed. These tests 
are described below. 

Leak checks . Leaks in the flow system were detected 
by pressurizing the system to 100 psig with helium and 
coating all fittings with a soap solution. The system was 
considered leak tight when no observable change in the 
Heise pressure gages could be detected over a span of 
several hours. Leaks in the manometers were detected by 
isolating 100 psig on each leg of the manometer and 
observing if any change in the fluid column occurred with 
time. The micromanometer was checked at approximately 
15 psig due to pressure limitations of the instrument. 

This process was repeated several times during the course 
of the investigation to insure system integrity. Pressure 
level manometers were only checked once (at the beginning 
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of * le testing) since the atmospheric side had to he 
connected into the flow system in order to provide a means 
of subjecting both sides of the fluid column to an equal 
pressure . 

A continuous check of the vacuum system pressure 
was provided by the vacuum gages present in the system. 

Differential pressure cell calibration . The 
Foxboro integral orifice, differential pressure cell was 
used as a secondary check of the primary flow measurement 
instruments. The 0.03^ and 0.0595 inch Foxboro orifices 
were calibrated by measuring the discharge of nitrogen 
from a cylinder of known volume. The procedure employed 
was essentially the same as that used by Magee (5^) except 
the cylinder, fabricated from a standard compressed gas 
bottle, was not immersed in an isothermal bath. A long 
pressure tap and a Chromel-Alumel thermocouple were placed 
inside the cylinder through a cylindrical steel plug 
located at the top of the bottle. The time history of 
conditions in the cylinder, of conditions at the outlet of 
the orifice, and of the D.P. cell output during the blow- 
down allowed the orifice constant to be determined by 
numerical integration. An empirical expression for the 
net orifice expansion factor, obtained from the ASME study 
of fluid meters (55) » was included in the data reduction. 
The results differed from the factory supplied calibration 
curves by a maximum of k per cent, for the smaller orifice. 
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The orifice constant for the 0.0995 orifice, used in the 
two highest Reynolds number runs, was determined from the 
factory calibration; limited capacity of the cylinder 
prohibited blow down runs of any reasonable time duration 
for this orifice. 

Thermocouple conduction error . A thermocouple 
attached to a heated surface changes the heat loss at the 
point of attachment, due to conduction along the wires. 
Local surface temperature is thus reduced. To provide a 
correction for this effect, the thermocouples were cali- 
brated in place for a number of different power levels 
without fluid flow with a procedure comparable to the 
preliminary study of Hess (48). A ceramic thermocouple 
probe was inserted through the exit of the test section. 

The conduction error was estimated by obtaining the inside 
tube wall temperature at several locations, far enough 
removed from the wall thermocouple to be uneffected by the 
thermocouple attachment. The wall temperature that would 
have existed at the wall thermocouple location if it were 
not present was predicted from a curve fit'*' of the probe 
temperatures. Correction of this value for the temperature 
drop through the wall yielded an estimated "undisturbed , H 
outside wall temperature which was then compared with the 
value measured by the thermocouple welded there. The 

1. Procedure FIT, reference 49. 
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ceramic probe was located by means of an external scale and 
reference positions marked on the probe. The zero refer- 
ence was obtained by null balance of a DC voltage signal, 
between the probe tip and a wall thermocouple located in 
the unheated region above the upper electrode. An estimate 
of the thermal expansion of the tube was included in the 
calculation of position. The maximum error in the location 
due to the difference in the thermal expansion of the two 
different materials was estimated to be less than 0.048 
inches at 800 °F . The correction obtained was expressed as 
a linear function of the wall temperature. 

Heat loss calibration . Data for the heat loss 
calibration were collected simultaneously with the thermo- 
couple conduction calibration data. In addition to the 
probe and test section wall temperatures, the current, 
voltage, vacuum pressure, and the temperature distribution 
along the glass vacuum vessel were recorded with each probe 
reading . 

The procedure for evaluating the heat loss was 
developed and was programmed by A. F. Deardorff. Local 
radiation heat loss was determined by subtracting the local 
energy loss due to axial conduction from the average 
resistive heating. The points were referenced to \/k inch 
intervals on the unexpanded tube and were corrected for 
expansion when heated. Small cylindrical, isothermal 
elements were considered to be located at these points. 
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Each dement was assumed to be surrounded by a 6 inch 
diameter cylinder centered on the point and divided into 
seven equal isothermal elements* A radiation network 
describing diffuse energy radiation from each inner element 
i was constructed and solved by electrical analogy (56). 

The element i was considered as a small gray body in a 
large enclosure * The resulting expression for the radiated 
heat flux, 



e . (T . )(T^ 

X xw xw 


8 

- I 


n=l 


L 

r. x ) 

xn n 


(E-l) 


was then used to evaluate the emissivity. For several 
different power settings, the emissivity data at each 
location ■, i, were correlated with a straight line least 
squares fit. 

Geometric view factors, F. .from the inner element 
to the seven outer isothermal glass elements, were 
evaluated from the defining expression for view factors 
(56). A numerical solution generated and programmed by 
Mills ( 57 ) » for evaluation of view factors in a non- 
isothermal annulus, was used to obtain these values. The 
angles in the definition for the view factor were evaluated 
using vector algebra. The view factors for the two end 
disks, at assumed constant temperatures of 530®F(, were 
evaluated using flux algebra, 


F iU = J 



F . 

m 


.126 
(E 2) 


An average temperature for each 5 inch isothermal glass 
element, T^, was evaluated from an interpolation of the 
measured glass temperature profile. The temperature at the 
point i on the test section was evaluated using the cor- 
rected temperatures and Procedure FIT (**9)» 

Adiabatic friction factor checks . Over fifty 
adiabatic friction factor runs were made, using air and 
helium, prior to the collection of the heat transfer data. 
The pressure taps were located 28.893 inches apart and the 
flow covered the bulk Reynolds number range from 750 to 
38,000. Comparison of the laminar results with the laminar 
prediction and the turbulent results with the prediction of 

Drew, Koo and McAdams (1) provided a check of the system 

♦ 

integrity and correct operation of the flow and pressure 
drop measurement equipment. Results of some of the later 
turbulent runs have been included in Figure 1. 


Heat Transfer Runs 

The heat transfer runs consisted of measurements of 
the test section current, voltage, overall pressure drop, 
pressure at the first pressure tap, wall temperatures , 
the mass flow rate, inlet and outlet bulk temperatures, 
the vacuum system pressure, and the temperature distribu- 
tion of the vacuum chamber wall. The mass flow rate 
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measurements consisted of the primary and secondary flow 
meter readings and the gas temperature and pressure down- 
stream from the meters. 

The successive procedure follows. Zeroes of all 
manometers and Heise gages were set under atmospheric 
conditions. The system was then pressurized to approxi- 
mately 100 psia. The fluid levels in the pressure drop 
manometers were alternately raised and lowered, by opening 
and closing bypass valves and vents, until the zero was 
consistently reproduced. This step was necessary to assure 
that only test gas was present in all of the manometer 
lines and in both legs of the manometers. The bypass 
valves were then closed and the system leak checked for a 
period of one hour. 

The Foxboro Differential Pressure Cell zero was 
then set, and the ice reference and zeroes of 4t »e manometers 
and meters checked. Correct operation of the Electronic -18 
over the expected temperature range was checked by means of 
a calibration feature of the ARU unit. 

The desired inlet flow rate was then set and an 
adiabatic friction factor run made to determine the 
validity of the pressure drop and flow measurements for the 
heating runs at these settings. The test section power 
leads were then connected to the power supply and the 
current adjusted to give desired wall temperatures. Care 
was taken to assure that the current for the setting was 
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near full scale on the ammeter. When steady state was 
obtained, as indicated by no further changes in the 
readings of the last three or four thermocouples located 
on the heated portion of the test section, the data 
described earlier were recorded in the data book. The 
meters used and the current transformer setting were also 
recorded. Overall test section resistance was computed 
from the voltage and current measurements, and a field plot 
of wall temperature prepared. After checking any readings 
that appeared out of line, a new power setting was set, if 
adequate supply gas remained, and the process was repeated. 
At the completion of the heating runs, the power supply was 
disconnected and an adiabatic friction factor run taken 
when steady state was reached. 

Data Reduction 

Except for changes in the calculation of the local 
energy generation, heat loss, gas stagnation enthalpy, and 
the mass flow rate, as necessitated by changes in the 
experimental loop, the data reduction formulae developed by 
McEligot (^9) were employed throughout the heat transfer 
study . 

Gas Bulk Stagnation Enthalpy 

The bulk stagnation enthalpy at a thermocouple, n, 
was obtained by an energy balance on a control volume 
comprised of the tube and fluid, as shown in Figure E-2. 



129 





Figure E-2. Control Volume for Calculation of Bulk 

Stagnation Enthalpy and Typical Axial 
Profiles 
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Trapezoidal integration was used to find the net difference 
between the energy generation and the radiation heat load) 
in the energy balance, 


H 

o ,n 




m-1 r 

+ q cond,n + ^ q g 



j + 1 


+ 





(E-3) 


where q ' . and q' . represent the energy generation per foot 

& y J r » J 

and the radiation loss per foot at the location j . The 
assumption was made that the small amount of energy 
conducted down the tube below the lower electrode all 
entered the gas. Calculations, based on temperature 
measurements at the electrode and on the wall below the 
electrode, showed that this assumption introduces a maximum 
error of less than 0.1 per cent in the bulk temper«ature at 
the lower electrode, for the air run at a Reynolds number 
of 3860. 


Conduction Heat Loss 

The heat loss at the lower electrode, q_ , was 

’ ^le 

estimated by considering the electrode as a thin circular 
radiating fin. The resulting expression (48) 


L le 


■vr 


le 


ji 2 0-k 


le e le 


(T 


le 


5T T 

? leoo 


^ + 4t ^ ) 
' * * l oo J 


1/2 


(E-4) 



was found to be closely approximated (for the present 
stainless steel electrode) by 
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q le = 0.0224(T le - T^) 90 < T < kOO'F (E-5) 


The latter expression was employed for the calculation of 
q. in the energy balance. A 100 per cent error in this 
approximation would cause a maximum error of less than 2 per 
cent in the bulk temperature at the lower electrode. This 
error would occur for the air data at a Reynolds .umber of 
3860 . 

The conduction heat loss was evaluated from the 
expression 


q c ond 
n 


-k(T n )Aj. 


dT 

n 

dx 


Procedure FIT 


(k9) 


was used to 


determine 


dT 

n 

dx 


(E-6) 


Local Energy Generation 

The local energy generation was calculated both at 
1/k inch intervals and at the thermocouple locations from 
the expression 

q* = I 2 R'(T.) (E-7) 

& » J J 

The positions were corrected for tube expansion. Procedure 
FIT (( 19 ) was used to obtain the temperature at the 1/k 
inch intervals from the corrected thermocouple temperatures. 
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The local test section resistance was evaluated as a 
function of temperature from the experimental expression 

0.12452 + 1.058 x 10“ 5 (T . - 80 *F ) 

R'(T .) « ; J S?/FT (E-8) 

exp 

The calibration for the variation of the resistance with 
temperature is described in Chapter VII. The major portion 
of the calibration , including a measurement of the axial 
variation of resistance for the test section, was performed 
by R. W. Shumway ( 58 ). The straight line correlation used 
for the resistance in this study is slightly lower, and 
approximately parallels the result obtained by Shumway ( 58 ). 
The new correlation was based on consideration of more 
recent data and the observed decrease of resistance with 
time, as indicated in Chapter VII. 


Radiation Heat Loss 

For the radiation calculations the test section and 
the vacuum chamber were considered to be a non isothermal 
annulus. The radiation heat loss at points located at 1/4 
inch intervals on the unexpanded tube, corrected for 
expansion when heated, was calculated from 


8 

W 0 " (T / - ^ F j,n T n 4) 


(E-9) 


where e.(T ) was evaluated from an expression of the form 
J J 
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e . (T . ) * a .T + b . 

J J J 3 .1 


(E-10) 


The a . and b . coefficients were obtained from the heat loss 
J J 

calibration. The calculation procedure is the same as that 
presented in the section--Heat Loss Calibration--with the 
exception that e^.(Tj) is now the known quantity , and "r',j ts 
the unknown. 

The radiation heat loss at a given thermocouple was 
ev.'*luated from a straight line interpolation of the values 
of ./T^ on either side of the thermocouple. 

Mass Mow Rate 


The mass flow rate was calculated from 

O 

m = sc fm 


A-VZ 


(E-ll ) 


P 


air density at 68°F, 1*1.7 psia (Brooks 
"standard conditions"), 

fluid density at calibration conditions, 
actual operating density at the flowmeter, 
scfm s reading from calibration charts. 


std 

P. 

p 


c 


Local Heat Flux 

The heat flux at each thermocouple was calculated 
from the expression 


<C,n 


q* - q* d^T 

= - & ! »_ - 4r : n - k(T )Aj. —a 

n * 2 


JID 


(E-12) 


exp 


dx 


i yk 


where l> is the expanded diameter of the test section* 
exp * 

o o 

Procedure FIT (k$) was used to evaluate d^T^/dx*". 

Constant Property Extrapolation 

The constant property extrapolation is explained in 
sufficient detail in Chapter V£1 and will not be repeated 
here. The procedure is illustrated in Figure which 
presents the extrapolation for the indicated thermocouples 
for the series of six heating runs taken at an inlet 
Reynolds number of 6800 . Uncertainty intervals, discussed 
in the following section, are also included* 

Uncertainty Analysis 

An uncertainty analysis was performed by the method 
of Kline and McClintock (30) to estimate the uncertainty in 
the experimental data* The uncertainties in the directly 
measured quantities are presented in Table E-3* These 
values are based on experience and manufacturers' specifi- 
cations, where available. 

The uncertainties of the extrapolated constant 
property Nusselt Numbers presented in Figures 6, 7, 8, and 
9 were obtained by plotting the uncertainty intervals of 
the data for the different heating runs on the logarithmic 
plots used for the extrapolation. From these values a 
straight line uncertainty band was estimated. The values 
of the band at a wall to bulk ratio of one were used for 
the uncertainty in the constant property Nusselt Number. 
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Table E-3* Heat Transfer Study: Uncertainty Intervals of 
Measured Quantities 


1. 1 cm 


Uncertainty 


Current 

Voltage : 

Fluke 
Wes ton 

Brooks Variable Area Flowme 
Model 1112A 
Model 150 

Foxboro D. P. Cell 
0.03^ orifice 
0.0595 
0.0995 

Wall Temperature 7 

Inlet Bulk Temperature j 
Thermocouple location 
Pressure tap location 
Resistance 
Diameter 


♦0.25# full scale 


+0.1% of reading 
To. 25% full scale 

or 

+0.0221 scfm c (.1% full scale) 

+0.004l4 scfm 

— c 

+0.5*1% of flow rate 
+_2 . 5% of flow rate 
£ 5 . 0 % full scale 

+2°F, 3/8% of reuding 
above 535*F 


+ 0.005 inches 

+ 0.002ft 

+0.0009 inches 


Pressure (Heise) +0.15 P»i 

Pressure Drop 

10*' Inclined Manometer +,0.01 in H 2 O 

60 M Vertical Manometer +.0.5 in H 2 O 
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Representative uncertainty values for two of the 
heating runs at a Reynolds number of 6800 are presented in 
Table E-4. Values are also presented in Figure 5* which 
illustrate the constant property extrapolation. 
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Table E-4k. Representative Uncertainty Values Tor Nusselt 
Number Data With Heating 


Re - 6800 , Helium 

Run 53 35 

(T /T. ) 1.090 1.379 

w b max _ — 
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APPENDIX F 

DERIVATION OF MEAT TRANSFER ANALYSIS 

The following derivation closely follows the method 
presented by Sparrow, Hallman and Siegel (35) and sum- 
marized by Kays (2). The major differences have been 
discussed in Chapter VI . 

The energy equation for turbulent flow in a tube is 
derived for the conditions: 

1. Constant fluid properties, 

2. Fully developed velocity profiles, 

3* Negligible viscous dissipation and axial conduc- 
tion, 

k. Steady, incompressible flow, 

5* No internal energy generation* 

The energy equation reduces to 

4: (i - y)fi(y) 2= = u + <7) (r-i) 

y^i - y) ay ay 3 * 

where 

p(y) = y^(y) + p 7 and ©<x + »y) = ^7 - 7^ (F-2) 

For a constant wall heat flux and an initial uniform fluid 
temperature, T^ , the boundary conditions are 
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^-“(x*,0) = -1 at the wall, - 3 ) 

By 

— ^(x , l) = O at the centerline, 

By 

and 

0(0,y) * 0 at the start of heating (F-^l) 

Use of Reynolds analogy, = £* , allows evaluation 

of |3(y) from the momentum eddy diffusivity distribution, 
equations 3-5 and 3-10. In terms of the present dimension- 


less variables 

this yields 


P(y> = \ 

tr — 

>[ y <L y • 

+ ycyl n r _ ?' 

y^tanh — 2 - y 1 + 2(1 - y) 

♦ k lF -5> 


y ! JL J L J 


where 





y* = 11 + 9«lll6 exp(-0. 272^9 Re x 10~ 3 ) 

+ [l5*83 exp ( ~0 . 9^98 Re x 10~ 3 )] 
and it = 0 .^ 225 * 

The last term in the expression for y* is negligible for 
Reynolds numbers above *1000. 

The velocity is determined directly from the 
defining equation for £ , 



(F-6) 


1*10 


Since equation F-2 is linear, 0 can be separated 

into the sum of a solution for the fully developed region, 

0 , and the entrance region contribution, 9 , 

OO e 

9-9 +9 (F-7) 

e oo 

For large values of x + , 9 approaches zero. This separation 

© 

of the dimensionless temperature is employed to obtain 
homogeneous boundary conditions fpr the entrance solution, 
as shown later. Caution must be used in the interpretation 
of the fully developed solution . The term "fully 

developed temperature profile" refers to a profile 
generalized in terms of variables which make it invariant 
with axial position. Such a profile is normally used to 
obtain the fully developed heat transfer solution (2). 

While 9^ describes the temperature for fully developed 
flow, it is not invariant with axial position and, thus, is 
not a fully developed temperature profile in the above 
sense . 


Fully Developed Solution 

At large x* , 9 = The energy equation becomes 


- y) a yL 


39 

(1 - y>P(y) -22 
3y J 


39 

* u + (7> -T 

ax 


(F-8) 


with the boundary conditions 


29 39 

- 3 s (x + ,0) - -1, —22 (x + ,l) * 

2 y By 


o 


(F-9) 


The requirements for fully established thermal conditions 


are 

or 


that the dimensionless temperature. (T 

w 

(T - T)/(T - T 

w w 


), and the heat transfer 


- T)/(T w - T b ) 
coefficient h 


» 


h 



(F -10 ) 


be invariant with the tube length (2). For a constant wall 
heat flux (qJJ,) these two requirements yield, in terms of 
the present variables, 




where 



An energy balance on the fluid shows 


(F-ll ) 


(F-12) 



k 

RePr 


constant . 


(F-13) 


Combination of equations F-ll and F-13 leads to 



RePr 


(F-l4) 


The dimensionless temperature, €^o, could be obtained 
directly at this point from equation F-8, after substitu- 
tion of F-l4, by integration. However, the alternate 
approach of a further separation of (^will be employed. 


This approach will yield a more tractable form of 0^ for 
use in the entrance solution. 

The general solution of equation F-lk is 

e oo= + F( *> (F - 15) 

where F(y) is an unknown function of y to be determined. 
Integration of equation F-13 yields 

0 b ~ RePr X (F-l6) 


where the constant of integration is zero by virtue of the 

boundary condition 9. I =0. Combination with equation 

b |x + =0 


F-15 for 9^ reveals that 


T - T, 


F< y ) - e oo- °b ■ 


(F-17) 




w 


l wk 


d9 

Substitution of the equations for ■■■■ ^ and 9^ into the 

ax 

energy equation (F-8) yields a second order, ordinary 
differential equation for F(y), 


- y) d yL 


(1 - y)^(y) ~ = u + (y) -j~ pr 

dyj 


<F-l8) 


with boundary conditions 


dF ( 0 ) , dF ( 1 ) 

_ - » _ 
dy dy 


= 0 


(F-19) 


Integration from 0 to y leads to 


1**3 


(1 - y)P(y) = -p(o) + f u + (£)(l ( F -20 ) 

•ly ' i 


The value of p(0) is obtained from the above equation with 

7 = 

P<°> “ / u+ <£ ,(1 -£><£ 

0 

The numerical result of this integration was used for £(0), 
rather than 1/Pr (see equation ,F-5) • The use of the former 
assured that the right-hand side of equation F-20 would 
approach zero as y approached one. If this criterion is 
not fulfilled, errors can result in F(y) due to a non-zero 
value of the right-hand side of F-20. The integral 
expression for 0(0) will reduce to 1/Pr only if the 
numerical evaluation of the integral of the velocity, in 
F-21, identically equals the correct value of the bulk 
velocity. 


A second integration of equation F-20 yields 


F(7) - F (0) = J 

0 


7rr 


jmp 


kyt r J 

+ RePr J u 


0 


Of 


RePr 


(F-22 ) 


Ikk 


where it is observed, from equation F-17, that 

-(T - T) 

Ky) - F(0) = > A ** 9 (F-23) 

\ vr 


Thus, an axially invariant form of the fully developed 
temperature profile is given by 


T 

w 


T 

w 



F(0) - F (y ) 

fToT - fllT 


(¥-2k) 


Equation F-17 shows that F( 0 ) can be obtained directly from 
the definition of bulk temperature, 


F (0 ) 




y)dy 

(F-25) 


The local Nusselt number is defined as 


Nu = 


hD 


(F- 26 ) 


T - T, 


Substitution of F -10 for h and F-25 for " y 11 yields 

q w r w /K 


Nu = 
OO 


fToT 


tF-27) 


Thermal Entrance Region 

The entrance region contribution is defined by 


®e * 0 - e co 


(F- 28 ) 



lk$ 

By the linearity of the energy equation, © also satisfies 

0 

equation F-l, 

_ a r a© I A a© 

- T — - — r~ (l - y)p(y) - u + (7) — 2. (F -29) 

y^(l - y) By L 3y J 3x 


By virtue of equations 

F-3 axid F-9, the boundary 

conditions 

are 



B0 a 
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— — (x , 0 ) 
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Substitution of the expression for €^, F-15» into F -28 
yields the initial condition 

© (O,^) = -F(y) (F-31) 

A product solution is assumed 

© e = R(7 )X(x + ) (F-32) 

Separation of variables leads to 

)2 + 

X = Ce~ A x (F-33) 

and 

~ (l - y)£(y) — (y)l + ^yici - y)u + (y)R(y) = 0 (F-3^) 

dy [ dy J . « 

dR( 1 ) _ 

» = 
dy 


dR(0) 

dy 


o 


(F-35) 


1^6 


The sep. ration constant is chosen so th.it 0 ^ will approach 
zero for large x . 

Equation F-3^ ami its boundary conditions form a 
Sturm-Liouville problem: are the discrete eigenvalues 

(though infinite in number) for which a solution exists , 
and yj^(l - y)u + (y) is the weighting function. The solution 
for R is expressible in terms of an infinite series of 
functions, R^, which are orthogonal with respect to the 
weighting function. R^ and are the ntf* unique combina- 
tion of a function and a constant multiplier which will 
satisfy the differential equation and boundary conditions. 

Thus, the solution for © e is given by 


\ 2 x + 

0 = V C R (y ) e 

e Li n n 3 


(F- 36 ) 


Application of the initial condition, F-31i at x = 0 
yields 


-F(y) = £ C n R n (y) 


(F-37) 


The orthogonality of R^ with respect to the weighting 
function, y^(l - y)u (y), allows evaluation of , from 
equation F-37, by using the principle of orthogonality 


(60) , 
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n 



- y)u + (y) 


F(y)R n (y)dy 


- y)u + (y) 


R ^(y)<*y 


(F- 38 ) 


Ik7 


The dimensionless temperature at any position is given by 
the sum of 0^ and © e 
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k ♦ . 

RePr x + 


1 ( 7 ) 


2 «► 


0° -) “x 

♦ y C R (7)e An 

La. 11 J| J 
11 = 1 


(F-39) 


Substitution of equation F-l6 for the first term of F-39 
and y = 0 yields the dimensionless wall.--bulk temperature 
difference 


vr 


- e b - 


T - T, 
w 1 

q"r /k 

V w' 


= F (0 ) + 


OO 
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n = l 


C R (O ) e 
n n 


1 2 ♦ 

"A,i X 


(F-*lO) 


Heat Transfer Results 

Combination of the definitions of the heat transfer 
coefficient and Nusselt number (equations F-10 and F- 26 , 
respectively) yields 


q M D 

Nu( * + > - TT~ - x b ) k 


(f-*u) 


w 


In terms of the present variables 


Nu(x 4 ) = 2 


0 - ©, 
w b 


Substitution of equation F-kO for © - ©, yields 

w o 


Nu(x + ) = •• 


2 + 


00 -A X 

F (0 ) + E C R (0) n 


n=l 


n n 


e 


(F-43) 


The fully develop©*! Xusselt number is obtained from the 


above expression by allowing x to approach infinity. The 
result agrees with the fully developed value predicted 
earlier by equation F-27» 

Ku oo “ rfoT (F - 27 

The normalised entry length X r usselt number, used for 
prediction purposes, is obtained by division of V by 
F-27 


Nu(x + ) 


Nu. 


00 


00 

1 ♦ [ A e 

, n 
na. I 


'An 


2 + 
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(F-W 


where 
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g 


C R (0) 
n n 

fToT“ 


(F-45) 


The equations presented above were solved 
numerically on an IBM 7072. The eigenvalue solution and 
other important numerical details are presented in 
Appendix G. A listing of the FORTRAN program is included. 



APPENDIX G 


DETAILS OF NUMERICAL HEAT TRANSFER SOLUTION 


Fully Developed Solution 

The fully developed solution follows in a straight 

forward fashion from equations F-22 and F-23» which were 

evaluated numerically using Simpson's rule. Prior to the 

integration, values of p(y) and u*(y) were calculated for 

the velocity profile in use at fixed y intervals (Ay - 

0.00125) and stored in the computer memory. (3(0), the 
C m 1 

value of rjj- + at y = 0, was evaluated from the expres- 
sion 


P(0) * 


RePr 


j 

0 


(y)(l - y)dy 


(G-l) 


The reason for using this equation is discussed in Appendix 

F. 

A comparison of the fully developed Nusselt numbers 
with existing formulations is presented in Table G-l. 

Eigenvalue Solution 

For the purpose of solution, equation F'-3^ was 
reduced to two first order equations, 

R' = S (G-2) 
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Table G-l . Comparison of Fully Developed Predictions 
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and 


s- = (7^= - -a yyr s 


)R 


1 - y 


P(y) 


(G-3) 


with R* (o) = O, R'(l) ~ C, The primes indicate differ- 
entiation with respect to y. P*(y) is obtained by differ- 
entiation of equation F-5* 


2 


The trial and error solution for pairs of 


and 


R was initiated by evaluating the coefficients of S and R 
n 

and storing them in the computer memory. Starting with the 

2 

initial trial value of ^ =0, the above set was solved 

simultaneously as an initial value problem. A fourth order 

predictor-corrector-technique with a Runge-Kutta method for 

starting was used (Hildebrand (39)) equations 6.6.1 and 

6.16.12). The process was continually repeated, 

incrementing ^ until a change in the sign of R'(l) was 

observed. Since the value of R’(l) varies in a damped 

periodic manner as ^ increases, the change in sign 

2 

represents the bracketing of a value of ^ for which 

R'(l) * 0. Linear interpolation was then used to predict 

the value of ^ for R'(l) * 0. The equations were then 

2 

solved with this value of ^ to find the corresponding 

value of R ’ ( 1 ) predicted by the differential equation. 

Linear interpol 1 1 on Wfcis again employed} with the two most 

2 

recent values of ^ and R'(l), and the process was continu- 
ally repeated to converge to the eigenvalue. Value e of 
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were then determined from equation F-38, using integration 

by Simpson's rule. The entire process was repeated until 

the first seven eigenvalues were found. 

Presented below is a listing of the FORTRAN program. 

It is divided into two main sections, the fully developed 

solution and the eigenvalue entrance solution. The choice 

is provided for a solution based on either the present 

velocity formulation or the velocity formulation used by 

Sparrow, Hallman and Siegel. Due to core limitations of 

the computer used, a separate program was employed for the 

McEligot, Ormand and Perkins solution. 

In the fully developed section, the eddy diffusivity 

and velocity values (for the profile chosen) are computed 

and stored for the desired Reynolds and Prandtl number. 

The fully developed temperature profile, (T - T)/(T - T c ) ; 

w w x 

F(y); and the fully developed Nusselt number, are then 

calculated in the order indicated, from the expressions in 
Appendix F. Simpson's Rule integration is employed in all 
integrations . 

In the entrance region, or eigenvalue solution 

section, the coefficients of equation G-3 are calculated and 

stored. Values of the velocity and eddy diffusivity 

calculated in the previous section are used where needed. 

The trial and error eigenvalue solution is then performed 

2 

to obtain the first seven values of ) , C , and A . The 

An’ n’ n 

entry Nusselt nuiiber is then calculated and normalized by 


153 

two different values, the fully developed Nusselt number 
and a Dittus-Boel ter formulation. The normalized values 
are then printed and new values of Reynolds and Prandlt 
number read. The program terminates if no new values for 
these two quantities are provided. 



n r * 








15 * 


pckthjn l i :• t I ::g of c^uii.^ program 


c 


O.Wvm »t P90GR 4 M l OiW RE Y.VJi. [)S W'lT' TURBULENT HE 
»»» AT TRANSFER 

0 I MENS I ON vP ( b-2 > « FPEY ( c ) OFK ( ow" 1 *vPlUG < -VC * ) WA 

* * <■ R ( V 0 ?\ 

t I MEN: 1 ON /. ( /•' 6 ) »Ff. < 20C- ) • £ < ?■:•(. ) ♦ <-( 2 . f ) ♦ & i 2 06 ) *S< EC 6 

*#* ),*.•'( VJf ) 

[: r M :_NS I ON f Y ( 1 -V 1 »JM I T ( 1 v. 4 ) • L.UMY ( 104 ) *A ( 9 ) *Gf I < 9 ) ,» 

* * “ M ' Jib .') 

0 :*-\;Noic\ Xr n ■:»< 3. h^(^i*oui* y p l o 5 ; : t> •:./'. > 

DIMENSION 0 P i ’. ( I 

DO W »\ i T R E A M , f- u 1. 1. V r M.VM';Pt!) FI.O W 0 0 1. 1 , j 7 ! 0 N 
4 READ £• RE ♦ 0C ♦ I EDDY 
3 FORMAT < 2F • 1 .’>.0 , I 1 ) 

CCspAriOOf FRiCT I^NFACTCK <,f- r I T ) /HL VJf IOf. FRICTION 
**< FAC TCP. 

cc should i-r :>[ r = to i. if Blaus us friction facto 

*** 0 IS USED 

LXPRt. S3 I ON rop FRICT H-N FACT oft A NJ T HcCRE T i C AL BU.il 
a** VELOCITY 

1 EDDY DLTCRM I N F 0 V E L O CITY PR 0 h I L i U S F. D 

i e: d o y * i mod i f i ;r o r l i c. >■ a r d t * 25 p a R p 0 w ?•: r a l • 

RC*D vDOCfpP 
40 )U FORMAT ( FI 0.0 ) 

FRIT--*. D791/<RF**.?5- ) 

FR1 t -CC^FR I T 

Y M A X * R E * S 0 P T r ( F o i r / ? , ) / ? . 

V5TH*SQRTF ( ? ./FR I T > 

PRINT :r. . • RE t YMAX 

300 FORMAT! IPX. lAHREYNOLDo HUMBER* * f- / .0 . 5X * iOUY PLUS MAX 

* * a * F 1 0 • 2 ) 

PRINT 3ul*CC 

301 FORMAT * 10X ♦41HFR I CT I F AC TOK/h’LAUS I S FRICTION FACT 

* ^ * po = , F 7 .4 ) 

NNN -NUiVSIER CF Y/.?y’ SPACES)*- IN MAGNITUDE TO DEL* USED 

IN INTEGRATION 

OF U+. NNN MUST* EVEN INTEGER MULTIPLE CF 80. THE 

RESULTING 

SPACING USED IN SOLUTION OF EIGEN VAL S PROBLEM A 

NO IN OBTAINING 

T - T W WILL BE 4DFL * AN u WILL* 8 DEL IN NUMERICAL INTO 
**» ORATION OF CN. 

NNN* SU'D 
A A “ N N N 
DEL* 1 • /AA 
J*NNN+ 1 
J E N 0 =• J / A +1 

Y R ( 1 J * 0 . S' 

INSERT FORM FOR FDD Y DTFFUSIVITY HE°E 


n r\ 
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GO TO ( 6 * SC #70 ) *r r :.'PV 
7 CONTINUE 

C vtcCIF i Ll. Rt l C 1 1 •*. fDT RKv-! II c 

6 OR I NT 1 

1 I ORMA r I iOX « 26HMCD I F I L ^ K L I C H*Rt> T PKGF I L l I 

(.. EXPRESSION FTP Yl 

h Y 

K t YN *K’M*C C I 

Yi = 1 1 *♦9.1 1 16/LXPf F < .4 7 ; u * u [. y ) 

IF ( RFvn- 7, ) I! # 1 ?♦ 1 ? 

1 1 V!..aVL i 1 S.BT/f :YP r .F ( *?44P*f<EYil) ) **4« 

1? CON T I Mt.'F 

C END OF EXPRESS ION FOR Yi. 

13 PRINT ] 4 • F K # Y L 

1 4 f ORMAT I 1 -X ♦ ?HK * *F 7 . 4 • VX • *H Y ♦ L* *F 7 . 3 1 
00 ICC I-i#J 
NM -fl 

YPLU5 ( I ) -Y R ( I >*YMAX 

If- { YPLUr,( I ) / Yl.-l G * )7fbtS 

7 fc*F_XPEFf v Pt.US ( I)/YL) 

T A NH - ( S -- 1 * / 6 ) / < P 1 • / b * 

GO TO 9 

8 TAMH-l.C 

? AA*RK»( YPLU«-( I >-YI*TANm) 

EDDY ( ! I = AA* (2.-YR ( I ) ) * ( 1.+2# * ( I 1* -YP tin * # 2 • > ) /6# 

D E R ( i » * I 1 ♦ -YR ( II) / ( 1 • + EDDY (!) } 
lOw YRINJ =YR< I) + DEL 

C I NT L : OR AT ION OF EDDY OlFFUSlViTY TO O'! TAIN VELOCITY 

YRDLL*LfcL»YMAX 
UPLOSI 1 ) = 0 • C 
ARE -0,0 
NIM-J-2 
K-? 

DO 2C0 I = 1 #NN * 2 

UPuOS ( IO -ARE* < DFR ( T ) + 4. *PL R < ' f 1 ) +15ER ( 1+2 ) >*YPPFL/3. 
ARE-UPL .UStK ) 

?00 <-K4l 

GC TO 229 

ENC OF MODIFIED RF I CHARD T PROFILE 
Sr ARROW tE T AL. EDDY DiF f-UoiVITY EXPRESSION 
SC D£Rli)=i. 

PRINT 2 

2 FORMAT ( 1 )X • 32HSPARRC /; ♦ HALLMAN# SlLGEL PROF ILF) 

Y PL 05 m *0. 

YR( 1 )~0, 

EDDY ( 1 I =0. 

YPOEL-DEL*YMAX 
PK- • 1 ? 4 * * ? # 

UPUJSJ 1 )*C 4 
UGtS=?#*YPDEL 
ARE-O. C 
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l-l 
u •- 1 

fr 6 N - ! ♦ 1 

v » | ♦ 2 

VP(N! -V9| I ) +''f u 

Vv' [ i S YR ( N ) 4 nr L 

YFLUSI \> «Y°LL'S < 1 I ♦ 

fp '„ ' { M i - vjjl:; j, ( M ; ♦ YfOci. 

7? l.v- i^'LU': ( L > 4-i UOt.S — w‘Pl. U *-» ‘x i)) / 2 • 

AA-PK*U« YPLI-r. ( N ) 

i; D D v ( N > * A A * ( 1 1 * /f.XPLF { AA ) ) 

D t K < N ) « * i • - Vf: t \ M / { • ♦EDti'Y IN)) 

A A-Rk^UOLS^YP^ JS ( »•« > 

fc‘DD v ‘ Mi *AA« ( 1*-1 ./rXPcF(AA) ) 

Dfc R t M ) - I 1 • - Y» TM ) / ( 1 • ♦ CODY ( M ) ) 

UPLUS ( <, ) = A '• F-f ( C'FR i I ) ♦ Nl+DCRCM) ) *YPDEL/3« 

D -U L 'l J.S< Y ) -» GF - 
IKI ABf>F (D) -.OC5 j if,* 1 i *99 
9<? ir»p) i c • i s • i-i 

] • GGE u ; ( J F LL S ( K ) ■ '.JOE 5» ) / 1 • *UGLG 
GO r O 77 

1 S I F ( YPI.US ( M ) -26 • ) 16* 1 ? *i / 

1 6 UG£.S = U°U .*5 I K > ♦UPL US ( < * -gPLuS ( L ) 

AR£-UPLUS( K ) 

K - r. ♦ 1 

L -L + l 
1 *1 f2 
GO TO A 6 
.17 88* • 36* Y M A y 
MH-sM+l 

YR(NN) - YR ( ,v ) 4 Of. l 

YPLUS ( NN1 »YPLOS f v 1 ♦YPOtlL 

DO 20 l. - 1 'J N • J 

f= DO Y ( u < ~b6* VP (L|»!l."iHO) 

OL R ( L 1 * i 1 » - Yfc < L » t / i 1 •♦i.oDV < L J 1 
N«L+1 

YPL5JS ( N ) * YPLUS ( l ) +YPDC.L 
2.y YR( M) *YR<L )^-DEl 
NN* J-2 

ARE*UPLUS(K) 

1 -K-fl 

DO '30 L*MtflN*2 

UPL'JS ( I ) =AR£4 ( DCR ( L I R ( Lt- i ) +u£R { !. +2 ) ) *?PtJL L/ j • 

AK£*UPLUS< 1 ) 

30 I -1*1 

GO TO 229 

C FJiNiD CF SPARROW * L T AL. t.DDY b I Ff US i V I T Y EXPRESSION 

C CALCULATION CF F(Y»-FU).CALlLT TF ( I J 

* 229 DfL2~2**DEL 

l- 0.0 


VA * K » s l « P I . U f > ( K > *( 1,-2 > 

) A /-:2K'<t.LZ 

CONOT**, *YUAX/ < r it*r>$ ) 

Vi V ; - I * 0 

B“L22=Z.«D P L2 
A A — * 4 / P R 
7. “0.0 
I = 2 

L - J/2-1 

DC 740 N = 1 *(..*? 

W f 1 I * *a>' + ( VA r v ( N ) *4 # *l'Ari ( ,\i-i ) ■♦ V A F’ ( N + 2 M * D E L 2 / 7 , 
w w - * ( l ) 

4.- I * ! + i 

r.0- Wf Jr' MD 5 

✓ 03*2 . »t‘0/Vf/rn 
01 PRINT 9. ;*VDP 

OJ F OR.VAK 1 JX#24HV»UI.K/VUULK. T H L'OUfc* T I C A L = • F 1 0 • ) 

2 * 0.0 

v»(i)=r.o 

L=JLNO-l 

DO ?41 1^1 ,L 

W( J)*-l.+W( ! )/BC 

*0.4* I ~ 3 

r . r:(i }*m n /( a •-?. >*( fidoy <►. )*rr + i.u 

4 1 Z = Z+OLL72 

£U JLNDI -O.w 

1*2 

L * JLND-2 
DO 250 N*]«L*2 

f Y U) *BT ♦ (EE { NM-4 **l UN + 1 ) ft.E ( N + 2 ) ) *OEL72/3. 
bT*!‘Y(M 

^ v 1*1+1 

r y { i ) 

CALCULATION 0^ TW -T/ TW-TC *CALLED D** ! T ( I) * F < C> * CALL 
*** ro T FO * AND NU- 

7-0*0 

Z2*2«»DELZ7 
I END* J/S+l 
DO 26 v 1 " 1 * I END 
DMlTt i )=rv( I )/FY( 1 1 N 0 * 

K*4*I~3 

DUMY ( I i *— F Y ( 1 )*<1 >*UPLUMK> 

60 2=7+77 
t . * J / 8 - 1 
T r 0*0*0 

DO 2 70 N * l * L » 2 

70 T f*C = TF 0+ ( DU^Y ( M ) +4 * *DuMY ( N + 1 )+DUMY(N+2) ) *72/3. 
TFO--TFO/BO 



r\ r \ 
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c 


c t- 1 vUl a . ion : n v / , i 'Mii no* : i t ,..al to r < y » 

0 0 2 h 0 ! = i * I f N 0 

; a ■ r / (; ^ryiiuoc' 

OUTPUT FO*UM:> FULL Y DEVELOPED FLOW 

up I M r »ri;NUt !T- 

2"- row AT ( 1 K.11H.V: Fi,'.l 1 TV L V[.i. opro- ,f 1 % ’ # *>x* f jHr ( 2 I •- • 

*** Pc 12. 4/) 

PRINT 291 

2 9 1 FORMAT ( 1 3X *4HV/P'/‘ ♦ 6 X * i ~;i T >v~ T / j w - TC • \> X ♦7r?F < Y / >• ) ♦ SX * 

* # * ?HY-» *7X *2HU-*># 3X 

1 • 1AHE0DY DIF,/ VIS./ j 
A*NNN/8'. 

DO 29 3 **1 • I END 

1 r.a#M»7 
L * 4 # M 3 

PRINT 292, YR< J ) ♦flMIT ( 2) , F Y ( Y ) ,YPLU3< I ) fUPLOS (L) * FDD 
*** YC) 

? 9 P F 0 R‘1 A T ( 1 3 X • J P F 5 . 3 * F J 2 • l . S X , 1 P r ] p . 4 f 0 P f 1 2 • 2 • r 1 0 • 2 • F 1 

*** T • 3 ) 

2/3 CONTINUE 

PNC Cir FULLY DCVlI.OPc.i- FLOW SOLUTION-START OF 
SOLUTION OF EIGFN VAL^t PRO&LEV- 
CALfULAT ION OF F ( M ) AND G(\) AND DERIVATIVE OF EDDY 

( * *** UI^FUCIVITY 

CO T0(299f3o0u»30Jl ) tU.DDY 
3201 CONTINUE 

C V C OIF UO RE l C H A R D T P K uF i L E 

2 99 t * 2 

V « ymax/yl 

2 “DELZ 

DO 400 **2,6*2 
B-Fxr- FF( V»7 ) 

D E N *- 6 + i , / B 
SECH-4./ ( DLN**2 , ) 

T ANh* ( <3-1 . /G ) /DEN 

Y Y ~ 2>*2, 

Y y y = ; » «• 3 # 

BDYP=(YMAX-YMAX»SECHJ*'(0.-li.«2 < fc##YY-2.*YYY) 

EOYP- ( L;.5YP + ( YMAX*Z-YL*TANtt) * ! - 1 1 • 4 i6 . -6 . *YY ) )^F'K/ 

# X it . 

O • 

OENO- < EDDY ? y ) 4-1 , /cjr ) 

E S i L ) *- 1 ./ { 1 .-7 ) ♦FPvp/DFNO 
6311 } - YMAX»!,PlUO ( K ) /DLNC 
L »L + 1 

402 2«2+DELZ2 
N’=P. 

Z* DELZ 2 
L * J/2+ 1 
004 1C K = 3 , L # 2 
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! F < V* 2-1 . « J 40 i *4 ■ > « 

4.1 L’*E>fTF(V*7 > 

L'r^.s'-if U/ii 
Sf’CHa^ • / < Df N»*3* ) 


> ; 
v S> 


T A iV ■< 

- ( r - 

1 

LO T 

0 4 v 


G lCH 



T A Nf i 

s i • : 


Y B 7 

*•2. 


Y / = 

Z 4*3 

• 


E *■• Y p ~ ( YMA> ~Y V A A *SF Of ) * i 6* - 1 1 *?» *‘VYY } 

LDY^- ( Lai p * ( YVA T AN'f) * ( - 1 1 . ♦ 16 . »7 -{., * Y Y i ) * «k / 

* ■» 4 

b • 

V=2*K-i 

rr NP = cTr ,v j m j # / np 

£ ( ) s £ l YP / c f; no 
G ( i\ ) c Y M A X * ( ‘PLUS < K » /U<~ NC 
2 *Z+uE 1. 2 i 

4 l w N=N41 

GO TO A 12 

tNO or mo:; into <u i chmk,, t pkgf i ut 

C ALULA 7 1 ON OF Li!'*) ANu GIN) l$!,Vo OOOSCa'.I T Al. fJX 

*** PRES SIGNS 

3wUw AA«Z6«/YM4X 

'•'•u.- ; . i ?/*#*? # ) * Y f * a x 
L~2 

7= DEL 7 

00 340 K. ~ 2 f fi * i? 


OENO”EDDY< V) 4i # /pr 
1F<Z-AA)3K #31U*32v 
31v d6Ba3B*UPLUS(K »»7 

tOYP- ( < YMAX*Z« ( 1 # -Z M/DENO 4uPLuE f *f ) ) *P,B 
.‘.OVPsFDYPa (!#+(-)* 4-p.Htj ) / tXP’i J ; ( Bbh ) ) 

GO TO 33 G 

3 2: F0YP«.36*YMAX*< > 

3 3w if S (L ) --1 •/ ( i # -Z » >EOYr/OF.NO 
GS < L ) *YMAX*UPU1$ i .< ’ /Dt.NG 
l. = L4 1 

34 2 Z*Z4DCL 7.1 
N*2 
DFL 7 ?. 

L-J/2 ♦ ! 

DO 408 F=3»L*2 

0 1 NO- E DO Y ( M ) 4 1 . / P p 
IF (Z-AA)350»3 5'j # 4 0S 
l i? • b <3 8 = 3 B * U P LU 5 ( K J « 7 

tr YP's ( ( YMAX*Z* ( ] • ~Z M /'UENO +UPLUS ( F ) ) 

E D Y P - E D Y ? * ( 1 . + < - ] . 4 2 b U ) / E X f N£ F ( n fc B > > 

GO TO 407 


l6o 


c 


V. 

c 

c 


c 

c 


c 


C 

c 



'4 «i ' 


;• ?6* vs‘£ v ■> 1 1 .-7.*z > 

r tv)- i.uy» /dev 

0 ( J - VMAX*LPLUC ! I / Df. Vj 

1 7 ♦ >.2Z 

’.-Nf] 


so ro ai; 

K\0 0<* r '»V»f iiOw * f ? u. • 

z =on.zz 

DO 41) Ns? ♦ L 

rtu )*-)./( i.-zuriNi 

Ail 2 -7 4 Df: 1.2.2 
c ( i ) - - 1 
v* I 1 ) - U • ^ 

I NOl: X = 1 

D E L £ G = S T E P SIZE 3 Y WhlCH EIGEN VALUE 15 ADVANCED • 
f 01. “ACCURACY TO WHICH bOUNDAK Y CONDITION AT V/Klv«l 
#* * IS DESIRED 

L I G T • D 
c i: i e g * • ; 

t OL=.C^3j 1 
J t* I Cj •= 0 

CALCULATION OF COEFFICIENTS r OR PRc D I C T OR-C CRRF.CTOR 
«♦-* EQUATION' 

iVPAR e DELZZ /24 • 

1>2*.*WPAR 
X«!?Y«*WPAR 
v - 4 7 • * W PA R 
T *9«*W'PAR 
WW*19.*WPAR 
XX-5*»WPAR 

rhnge kutta start or solution 

4 ? (' f: ( 1 ) " D « C 

r m * u o 

SPl 1 )*-EICH-< D^RID 
Ot LH~CELZZ 
DO 43C I *1*3 
N = I+1 

0 a - 1) t L W * ( E ( l )*S( I } fCIo*C;( i )*Rl I) ) 

GG--CELH*U*MN)M5< I ) ♦ • E *0 ) ♦£! ( N ) * t R ( I I + .'j^DLLZ 

{ 1) ) > 

OOG* ( LS < N ) * ( S m 5 *CN n 

GGQ^-DElH* ( QQO-f f I C *6S ( N ) * < S( i >4,$*01LZ7.*S ( I )4-«25-»DL 
*## L Z Z * G J ) 

GOOO* ( E ( N ) * ( 5 ( I ) 4000 » > 

G0GQ S ~DELH4 (QGQO+E I G*G ( N »*((<( I )+DCLZZ*S< i ) + «5*f)f.LZZ 


» ** #00 ) ) 

R ( N )*R! i >40ELZZ*5 ( 3 1+U£LZZ*( J + GO+QCG ) /6* 

G < N ) ~ £ ( I > -f ( Of 2 •* ( OG*Gv*t. ) 4 GO GG) /6* 

A 3 r 5. P ( N ) « - f E < M * S I N ) 4E I G*U ( N ) *R < N n 

CONTINUATION OF SOLUTION bV MH ORDER PREDI C T OR-COR 
*## RECTOR 


I- * J ’ ft [) - l 
u'044 J 

P f 0 l TOP 

I? < ft i - R ( N- i ) * ft » ** ( * - ' )•*/*:( ft - » ) ♦ Y * M N - 3 ) - T * c , ( 
( t i * * s < ’ - ] > 4 * ■ ' p ( ft - I > - a * , » ( ft - 2 ) ♦ y * 7 , r ( f . - ? > - j 
' *M . ) *- ( i ( » « r. f \ h i i ( n » j 

CC.'PPE c TCP 

• v t N ) fc tt ( N - i ) ♦ 1 • G ( ft ) 4 >/ v, A w t ft- i ) • >x* -}{ lt-‘/ ) 4 f. | f,- 

• j ! N 1 * S { N- l ) 4 T * i ( \ ) 4 r. «. •> „■ ft t ft- / i +' 

* r « f'AP 

4 4 .. O P < ft ) * - ( L ( N ) » f> t ft j f f i 0 * <J ( ft . -4 i- 1 < h; ) ) 

N* JF \D 

^ { N ) - R ( l. ~ * ) ♦ * % { *•> - j ) - x * i, < ft - 2 ) 4- y * o ( N - 3 ) - T * S ( 
•' (\ y - ■: >' N- 1 ) ♦ ' « r » n • ft - l ; * I, ' • ( ft -v > -* Y •< r>P f - 2 > - T 


r. r' 


c 

c 


) r 0< b. j«4*:0.44.; »4Pc * • J ftftt x 

<»*-o rrtr»*rir 

n:sT*s( jl np j 

F. r^-C I G+l>f„ L F v « 
f ftCF.X - 3 
Gu Tc; 42 ) 

40 RAT 10«5>l JF. ftp > / 1 f. 5T 

!r "v A ’ I C) 4C > » .4b*« 

“to IF (AcSf ( S ( JF.riO M - Tm ) t>0 J * Jo » 490 
4V.J i'OG-UC 

(. A T S ( JEM)) 

T OP -• S ( JlftftD ) *»E l - T ! S T * t_ J G 
T n T M ( JFMP ) - r f ■ r r 

.g^a jl rAanr { f ig-k i r n > 

I t r i SMAl '.-•CC»',0C.G5 ) fn of »4Gi 
491 COST J N'JF 

1 F ( A33F f TH TM • -./Vft V'Jo4 > ;>1,0 *4 94 *494 
494 CONTINUE 

I F ( A t> : 1 1 i J r< ; • • sj 0 ( a 2 * l? o «. ♦ '♦ > v • 4 9 b 

4 9‘ •; I G ; - 1 OP / ‘['P' 

L tGP=DOG 
T i.ST=CAT 
i ftDF X = 4 

00 TO 42. 

CALCULAT ION OF ON, ft'; i f 0U V V l J J *D> , I T ( I > 

*** MEAN INC 

5» Of I«] 

JL l C> = J t I 0+ 1 

2 * 0.0 

» ' w 1 o *•* = 1 t J E ND « 2 
N«2*M- } 

0 1 J« R ( M I * U F'LU S ( N j * ! i • - 1 ) 

Dl'VY ( I ) =-F v ( I ) *;k v 

C)V I T { I hR( ,M , #r; UV 

2 - 2 -*Z? 

bl 0 T = T + 1 

JEN =0.0 


\rc 
#' • 


ft - 4 ) 

* J.P ( N 41 


U'n'AH 
P ( ft- 1 I •/ 


ft — 4 ) 

* '••TM n.-*, ; 


ftVf ft F V, 
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c 


I'TMO.U 
t. 'J/6- 1 

DO b2C 1 = 1 « L ♦ 2 
0 t N * Dfc N* ♦ ( O'JV V ( ! > 4 /« * » o ON : Y ( 
rt r y »i> T m ♦ ( dv it ( i i ♦ a , •* u. . j i 
r -c t \ / f> t ** 

A < jl !G i *c*p < l > / rr - 


♦1 5 4-CuMV ( ; *2 ) l*/.Z/3« 

: + : ; ♦u-mi t < !+?))*///)• 


of i ( jr i '• > -r ic 

PRINT 52 5 » f I G « o » s : ;f \D J • A ( jr. ; G ) 

5 2* FORMAT (1 X * 1 2>‘r I Of. t{ VALUE.* 3 '< ♦ 3HC N = • 1 PF 12.A 

* « * ,5>. , 3 HR * f 1 ) = ♦ PF 

1 l 2 • A ♦ 3 X * 3HA N = ♦ 'P r 17 .6 > 

P'Mf\ T 5 3‘- 

5 j - F JrtVA I I 1 3X *4h y , <• w « 3 X • 1 3HL" I Or N FUNCTION (R>) 

L =NNN/40 


DO 541 N r 1 * Jc. NO »L 
I » 3 

PRINT 540 ♦ YR ( 1 * • R I N ) 

5 4 : F OPM A T < i *» X » F 5 • '» ♦ 6 X • 1 P 1 1 2 • 4 ) 

541 CONTINUE 
I.NPFXa2 

f ! o- r ig^pfufo 

IF i JtIG-7l420*0 e ;»Dt55''i 

CALCULATION AN*NU/Nu FULLY Ot Vf.t.PPh D. 
5;.. X Y Z •• 0 * 

X =C 

L>B»FONU/ I •021* ( Rl *». 3 * * i P.'<*» *4 > ) 


DO 570 Ns.;, no 

K='>1 

X PLUS I N) *XY2+S, 

:> um *o*o 

DO 560 t. *1*7 

5 5 J SUM* SUM* .( L > /FXPEF ( Gf l ( L ) * XPUJSI N ) ) 

RMJ I N) * 1 . / U • + Si; M ) 

RDIMN) = RNL'(N) * D ? 

I F { RN'J ( N ) - j . 0U2 5 ) 5 a 5 * PH 5 ♦ 5 70 
5 70 XYZ*xr>u;£(Ni 
5 RS PRINT SBC 

SHJ FOkMAT ( ISXtOHX* ,3X*21T'NU/NU FULLY Ot VFLOPF.D ♦ 1 X * 7* ,Nu 
*** /NUDB) 



DO 591 N*1#K 



PRINT p 9; #yPL'JS IN) « 

ROu ( K ) , ROB ( N ) 

SVC 

FORMAT! .: Cx.cPFL.3 *F 

1 3 • 4 * F 1 V , 4 ) 

S 9 J 

CONTI NUt 



PRINT 592 


5 92 

FORMAT I 1H1 > 



GO TO 4 



END 



APPLMMX II 


TABULATED IILAT TRAiNSKER DATA 


Aoinonc 1 ature : 


T /T. 
w b 


Nu 

Re 

Pr 

+ 

q 


DB 


balk static temperature, # F 

wnil-to-bulk temperature ratio, °R/°R 

normalized parameter from Pittus-Boelter correla- 
tion, 0.021 Re 0 * 8 Pr0.4 

O 

Reynolds number , ~~ , properties based on local 


Prandtl number, 


c [l 
»\ * 


, properties based on local T. 


he£»t flux parameter, 


l" 

w 
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NOMEiVCLATURE 


English Symbols 

A quantity used in calculating Nussolt numbers from 

equation 6-5 

C coefficient in eigt,;' v.-»iue solution 

n 

C turbulence velocity correction? defined by equation 

A-6 


c 

P 

I) 

e 

F . 

in 

F(y) 


S 


c 


11 

o 

h 

k 

o 

m 

P 

A P 

q 

q' 

q" 


specific heat at constant pressure 
tube diameter 
emissivi ty 

radiation view factor from i to n 

function arising in fully developed solution, 

T - T t 
b 

q"r /k 
w w 

dimensional constant, e.g., 32 . 17^ ( 1 bm) ( f t ) / 

( lbf ) ( sec^) 

stagnation enthalpy 

q" 

heat transfer coefficient, h = ^ — ■ y" 

w b 


thermal conductivity 

mass flow rate 

pressure 

pressure drop 

heat transfer rate 

lineal heat transfer rate, 

heat flux 



l68 
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r 

It 

a 

It 

S 

1 

u 

a ' 

u* 


V 

w 


b 


III 


x 


y 

+ 

y i 

z 


tube radius 

eigenfunction used iu the separation of variables 
gas cons t . mt 

° • 
function defined b\ equation G-3 

temperature 

velocity in axial direction 

turbulent velocity fluctuation in the axial 
direction 

shear velocity, 

r 

bulk veloc ity , 

uncertainty interval of the quantity m 
axial distance from the start of heating 
radial distance from the tube wall 

Reynolds number dependent coefficient in correla- 
tion of eddy dil'fusivity for momentum 

compressibility factor 


/ 


u(r)r<ir or m/M Ur 




w 



Greek 

e<y) 

7 

A 

€ 


Symbols 

(Tm 1 

total thermal diffusivity, - — + 77 - 

]J Pr 

intermi tt ency factor 

effective center displacement, equation A-10 
eddy diffusivity; C for momentum; Cjj for heat 


0 


* 



T - 

dimensionless temperature, 7*i7JT 

constant of velocity profile 
eigenvalue 
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V 

p 

r 


absolute viscosit\ 
kinematic viscosity 
don si t> 
shear stress 


Dimensionless Groups 
f friction factor, 


T jr 
* w < 

P v ii 


iNu 


Dr 


Xus.se] t number, 7-— 

7 k 


Prandtl number, 


^c. 


Re 

+ 

u 


4 


heat flux parameter, *37 


11 

¥ 


ai 

A o 9 1 


P V o U 


Reynolds number,'—-*—— or ^m/TcD^l 

r 

, . , , u 

velocity parameter, 


b 

+ 


bulk velocity parameter, 


dimensionless axial distance, — 

7 r 

w 


dimensionless radial distance, 


w 


radial distance parameter, ^77— = ¥r 

V 2 r w v 


Subscripts 

b bulk, evaluated at bulk static temperature 


cond 


conduction 


roll") 
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s 

i 

w 

OO 

<L 


entrance 

generated 

inlet 

wall 

fu~ly developed 
centerline 
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